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the  project.  Following  the  Introduction,  in  Part  I,  the  research  undertaken 
at  Purdue  Univesity  on  analytical-computational  studies  is  presented  along 
with  relevant  publications.  In  Part  II,  the  research  undertaken  at  Imperial 
College,  London,  on  expci  imental  studies  conducted  to  establish  various 
aspects  of  the  influence  of  free  stream  turbulence  (isotropic  and 


nonisotropic)  on  low  speed,  fully  developed,  turbulent  boundary  layer  is 
presented  along  with  relevant  publications.  An  addendum  to  the  Final 
Report  includes  two  appendices,  one  on  the  large  eddy  interaction 
hypothesis  and  the  other  on  structural  similarity  between  velocity  and 


temperature  fields. 
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1.  INTRODUCTION 


The  general  problem  of  interaction  between  two  distinguishable  fluid 
flows  is  of  fundamental  interest  and  also  of  common  occurrence  in  nature 
and  engineering.  The  problem  becomes  considerably  complex  when  one  of 
the  flows  is  vortical  or  turbulent  even  when  both  flows  are  incompressible. 
The  presence  of  free  stream  turbulence  in  a  wall-bounded  flow  provides  an 
example  of  such  a  complex  flow.  It  can  be  expected  to  affect  both  skin 
friction  as  well  as  heat  transfer,  and  also  transition  in  the  boundary  layer 
when  otherwise  the  boundary  layer  may  have  remained  laminar.  The 
physical  processes  involved  in  such  interactions  are  not  fully  understood 
and,  therefore,  prediction  schemes  have  remained  generally  ad  hoc  even 
for  the  simplest  type  of  interactions. 

The  effects  of  free  stream  turbulence  are  of  interest  in  severa'  internal 
and  external  flows  in  engines  and  flight  vehicles.  Those  flows  involve 
additional  complications  such  as  pressure  gradient,  wall  geometry  and 
various  means  of  cooling,  for  instance  film  cooling.  Also,  there  can  be 
complexities  in  the  nature  of  free  stream  turbulence,  inhomogeneity  in  the 
large  scales  of  turbulence  being  a  common  occurrence  in  combustion 
devices  and  atmospheric  flows. 

A  body  of  literature  representing  the  current  state  of  knowledge  in  the 
subject  area  can  be  found  in  Refs.  1-12.  In  general  it  is  found  that  there 
has  been  no  systematic  study  of  the  effects  of  inhomogeneity  in  free  stream 
turbulence.  An  interesting  study  in  this  connection  is  reported  in  Ref.  13 
on  peaky  turbulence. 

In  the  study  of  interaction  of  turbulence  flowfields,  such  as  free  stream 
turbulence  (FST)  and  boundary  layer  turbulence  (BLT),  the  physical 
processes  may  be  fully  understood  only  by  an  appeal  to  a  combination  of 
the  classical  structure  as  depicted  by  statistical,  point-wise  information  and 
the  more  modern,  quasi-deterministic,  physical  structure.  Data  sets 
available  on  classical  structure  are  not  always  without  ambiguity  and  in  no 
case  include  pressure  fluctuations  in  the  flowfield.  There  are  no  flow 


visualization  or  quantitative  data  on  physical  structures  in  a  boundary  layer 
with  heat  transfer  or  free  stream  turbulence. 

It  may  be  worth  pointing  out  that  pressure  fluctuations  are  of  interest 
in  turbulent  flow  interactions  even  under  the  conditions  of  incompressible 
flow.  Unfortunately  there  are  no  data  sets  that  relate  free  stream  pressure 
fluctuations  with  even  wall  surface  pressure  fluctuations.  The  main 
difficulties  are  two-fold:  (i)  obtaining  pressure  fluctuation-data  with 
spectral  information  included  and  (ii)  extent  of  data  required  in  order  to 
relate  time-dependent  variation  of  pressure  to  (a)  the  velocity  scales  of 
interest  otherwise  and  (b)  the  dynamics  of  turbulence  production  adjacent 
to  the  boundary  wall  (Ref.  12). 

In  analytical-computational  studies,  understanding  in  terms  of  physical 
structures  and  their  dynamics  may  become  feasible  as  there  comes  about 
further  progress  in  large  eddy  simulation  or  solution  of  Navier-Stokes 
equations  (without  time-averaging)  and  the  needed  computational 
capability.  On  the  other  hand,,  one  can  try  to  obtain  the  best  out  of 
classical  description  of  turbulence.  A  significant  proposal  in  this  regard  is 
John  Lumley's  so-called  rational  description  of  turbulence  (Ref.  14,  15).  It 
is  interesting  that  turbulence  models  that  can  be  generated  from  it  have 
been  utilized  to  extract  physical  structural  details  in  various  contexts  (Refs. 
16-19).  Predictive  tools  may  also  be  developed  based  on  that  description  of 
turbulence.  This  is  the  basis  of  the  analytical  computational  work 
undertaken  in  the  current  project. 

1.1  Objectives  of  the  Investigation 

The  original  objectives  included  both  experimental  (at  Imperial 
College,  London)  and  analytical-computational  studies  (at  Purdue 
University). 

The  flowfield  considered  was  a  simple  one:  incompressible  flow  past  a 
flat  plate  under  conditions  of  heat  transfer,  with  no  external  pressure 
gradient,  when  the  boundary  layer  is  fully  turbulent  and  the  free  stream  is 
turbulent  and  homogeneous  or  with  a  few  dominant  large  scales  in  it. 


The  objectives  in  the  experimental  studies  were  to  obtain  measurements 
of  turbulence  quantities  and  their  structure  in  the  classical  sense.  In  the 
analytical-computational  studies,  it  was  proposed  to  adopt  the  classical, 
statistical  description  again.  However,  turbulence  was  sought  to  be 
modelled  on  the  basis  of  John  Lumley's  rational  description  of  turbulence 
and  its  adaptation  at  Purdue  University  in  the  form  of  what  has  come  to  be 
called  the  large  eddy  interaction  hypothesis  (Refs.  20-24).  The  modelling 
was  to  be  developed  for  obtaining  both  the  momentum  and  stress 
distributions  and  the  heat  transport  while  the  illustrative  calculations  did 
not  include  the  latter. 

1.2  Main  Accomplishments 

1.2.1.  Publications 

1.  S.N.B.  Murthy  and  S.K.  Hong:  Turbulence  Boundary  Layer  with  Free 
Stream  Turbulence,  AIAA  90-1503,  AIAA  21st  Fluid  Dynamics,  Plasma 
Dynamics  and  Lasers  Conference,  June  1990.  (Attached  at  the  end  of  Part 
I  of  this  Report). 

2.  V.  Bhaskaran,  O.E.  Abdellatif  and  P.  Bradshaw:  Effects  of  Free 
Stream  Turbulence  on  Turbulent  Boundary  Layers  with  Convective  Heat 
Transfer,  Seventh  Symposium  on  Turbulence  Shear  Flows,  Stanford 
University,  August  1989.  (Attached  at  the  end  of  Part  II  of  this  Report). 

3.  V.  Bhaskaran  and  P.  Bradshaw:  An  Experimental  Study  into  the 
Wake-Boundary  Layer  Transition,  10th  Australian  Fluid  Mechanics 
Conference,  December  1989.  (Attached  at  the  end  of  Part  II  of  this 
Report). 

4.  S.N.B.  Murthy:  Turbulence-Mean  Discontinuity  in  Supersonic  Flow, 
Meeting  on  the  Physics  of  Compressible  Turbulent  Mixing,  Princeton 
University,  Princeton  NJ,  October  1988.  (Attached  at  the  end  of  Part  I  of 
this  Report). 
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1.2.2.  Scientific  Advances 

1.  Application  of  large  eddy  interaction  hypothesis:  A  four-layer 
boundary  layer  model  has  been  developed,  along  with  a  numerical  code  for 
prediction  of  the  flowfield,  for  an  incompressible,  zero-pressure  gradient, 
fully  turbulent  boundary  layer  with  moderate  levels  and  scales  of  (grid¬ 
generated)  free  stream  turbulence.  The  Hancock-Bradshaw  measurements 
at  Imperial  College  have  been  reproduced  in  most  parts  utilizing  a 
skewness  and  a  damping  factor.  A  converged  solution  can  be  obtained  in 
400-600  iterations. 

2.  Logarithmic  region  of  boundary  layer:  The  numerical  work 
described  in  item  1  provides  the  basis  for  a  conclusion  that  a  boundary 
layer  type  solution  in  the  presence  of  free  stream  turbulence  requires 
necessarily  a  region  with  mean  flow  governed  by  the  so-called  log  law.  In 
the  current  predictions  it  has  been  essential  to  assume  the  existence  of  a 
layer  governed  by  log  law;  the  extent  of  that  region,  of  course,  is  a  part  of 
the  solution.  In  other  words,  the  presence  of  free  stream  turbulence 
permits  a  variety  of  mean  flow  profiles,  when  only  initial  and  boundary 
conditions  are  specified,  and  the  physical  flow  can  then  be  selected  only  on 
specifying  the  existence  of  a  region  governed  by  log  law. 

3.  Skin  friction  in  a  turbulent  boundary  layer:  In  the  experimental 
studies,  it  has  been  found  that  the  skin  friction  in  a  turbulent  boundary 
layer  on  a  flat  plate  decreases  in  the  presence  of  the  wake  of  an  upstream 
mounted  cylinder  representing  an  axial  turbomachine  blade,  in  contrast  to 
the  increase  in  Cf  produced  by  multiple-rod  grids.  Evidently  the  primary 
effect  is  the  reduction  in  "free  stream"  velocity  seen  by  the  boundary  layer, 
as  a  consequence  of  the  velocity  defect  in  the  cylinder  wake:  interaction 
between  the  turbulence  fields  of  the  wake  and  the  boundary  layer  are  less 
important.  The  details  of  the  flow,  particularly  in  the  present  case  where 
the  cylinder  axis  is  inclined  to  the  plane  of  the  boundary  layer,  are  of 
considerably  interest. 


4.  Skin  friction  coefficient  and  Stanton  number:  From 

experimental  studies,  it  is  found  that  skin  friction  coefficients  and  Stanton 
numbers  were  measured  in  a  flat  plate  boundary  layer  beneath  grid¬ 
generated  free  stream  turbulence,  with  a  maximum  rms  intensity  of  8 
percent  of  the  free  stream  velocity  and  a  maximum  length  scale  of  6  times 
the  local  boundary  layer  thickness.  The  increase  in  skin  friction 
coefficient  agrees  well  with  the  correlation  curve  based  on  the  variable. 
Vu?  /  Ue/(Le/8+2),  proposed  by  Hancock  &  Bradshaw  (1983).  Stanton 
number  increases  do  not  seem  to  follow  the  proposed  skin  friction 
correlation  or  any  simple  multiple  of  it,  indicating  that  the  Reynolds 
analogy  factor  is  changing  significantly  -  more  than  in  the  work  of  Blair 
(1983).  A  change  in  the  above  empirical  variable  to  improve  Stanton 
number  predictions  is  tentatively  suggested. 

5.  Influence  of  moving  rod  arrays:  The  wake  of  a  span  wise  rod 

mounted  ahead  of  the  test  plate  to  simulate  turbomachine  blade  wakes 
produced  a  reduction  in  skin  friction,  calling  into  question  the 
interpretation  of  published  results  from  moving  rod  arrays. 

1.3  Outline  of  Report 

The  report  is  divided  into  two  parts.  Part  I  presents  the  work 
undertaken  at  Purdue  University.  Section  2  deals  with  turbulent  boundary 
layer-free  stream  turbulence  interactions.  The  method  of  calculating  a 
flowfield  utilizing  LEIH  is  illustrated  in  the  attached  publications,  wherein 
the  flowfield  measured  by  Hancock  and  Bradshaw  (Refs.  25-26)  has  been 
chosen  for  prediction.  Section  3  then  deals  with  the  possibility  of  invoking 
classical  structural  similarity  between  turbulence  intensity  and  temperature 
variance  in  different  forms  across  a  boundary  layer  during  heat  transfer. 

The  conclusions  from  the  overall  effort  and  the  type  of  extensions  that 
appear  significant  are  discussed  in  Section  4. 


Part  II  deals  with  the  effort  accomplished  at  the  Imperial  College  of 
Science,  Technology  and  Medicine,  London,  U.K.,  under  a  subcontract  for 
experimental  studies  in  the  subject.  This  part  consists  of  the  report 
received  from  Professor  Peter  Bradshaw,  the  Principal  Investigator  for  the 
project  at  Imperial  College. 


PARTI 


SUMMARY  OF  PURDUE  UNIVERSITY  WORK 
ON  FREE  STREAM  TURBULENCE 
AND 

HEAT  TRANSFER 


S.N.B.  MURTHY 
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1.1.  LARGE  EDDY  INTERACTION  HYPOTHESIS 


The  general  background  to  the  development  of  the  hypothesis,  LEIII,  is 
provided  in  Appendix  I. 

A  number  of  features  of  the  hypothesis  and  its  application  to  boundary 
layer  problems  are  are  also  presented  therein. 
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1.2.  BOUNDARY  LAYER  PREDICTIONS  WITH  FREE  STREAM 

TURBULENCE 

The  asymptotic  matching  of  the  boundary  layer  regions  with  the  free 
stream  has  been  discussed  by  T.  Tennekes  (Ref.  27),  K.S.  Yagnik  (Ref.  28) 
and  G.L.  Mellor  (Ref.  29),  in  all  cases  in  the  absence  of  turbulence  in  the 
free  stream.  The  objective  has  been  to  establish  the  connections  between 
the  so-called  inner,  viscous  layer,  the  outer  layer  and  the  free  stream. 
Various  types  of  assumptions  are  intro^'ced  in  the  course  of  asymptotic 
analysis  and  matching  of  layers.  For  example,  a  key  assumption  of  the  case 
of  G.L.  Mellor  is  that  the  Reynolds  stress  is  of  the  same  order  of 
magnitude  in  the  inner  and  the  outer  layers. 

An  important  conclusion  of  G.L.  Mellor,  then,  is  that  the  outer 
Reynolds  stress  asymptote  of  viscous  region  (the  inner  layer)  must  behave 
just  the  same  as  the  inner  Reynolds  stress  asymptotic  of  the  defect  region 
(the  outer  layer).  In  fact,  therefore,  the  asymptotic  analysis  provides  a 
basis  for  a  region  governed  by  the  log-law,  which  was  rationalized  by  F. 
Clausen,  based  on  the  postulates  of  L.  Prandtl  and  von  Karman,  namely 
the  velocity  distribution  and  the  reduction  in  velocity  from  that  of  the  free 
stream  must  both  be  functions  of  wall  shear,  distance  from  the  surface  and 
fluid  kinematic  viscosity  in  the  first  case  and  overall  thickness  of  boundary 
layer  in  the  second. 

On  including  free  stream  turbulence,  several  considerations  arise  on 
amount  of  the  intensity  and  the  length  scale  of  such  turbulence.  Both 
intensity  and  length  scale  vary  along  the  flow  depending  upon  the  manner 
in  which  the  free  stream  turbulence  is  generated. 

Returning  to  the  problem  of  asymptotic  matching  the  free  stream 
velocities,  Reynolds  stresses  and  coordinates  may  again  be 
nondimensionalized  based  on  a  velocity  uu  ,  a  characteristic  turbulence 
velocity  ut  (which  may  be  the  same  as  ux  ,  the  friction  velocity)  and  a 
length  scale  i  . 


Two  small  parameters  that  appear  in  the  describing  equations  for  the 
boundary  region  are  the  following. 

£  =  and  e2  e  =  — —  , 
u0  u0 1 

where  v  is  the  kinematic  viscosity.  G.L.  Mellor  then  introduces  a  length 
scale  for  the  defect  region,  namely 

A/=  £  t  .  (2.1) 

It  may  also  be  noted  here  that  £  and  e  are  in  the  nature  of  ratios  of  (length 
scales)  in  the  viscous,  the  defect  and  the  free  stream  layers.  With  free 
stream  turbulence  of  a  given  scale  4  one  can  construct  a  length  scale  for 
the  free  stream  layer  that  becomes  coupled  to  the  boundary  layer,  namely 

Af  =  £*  4  -  (2.2) 

where  £*  is  given  by 

o*  =  Uf 

u,  (2.3) 

with  Uf  =  1/3  •  Ufy  ,  the  representative  turbulence  velocity  of  the  free 
stream. 

It  is  then  necessary  to  distinguish  between  the  following  cases: 

(i)  Af  «A,, 

(ii)  Af  =  A/ ,  and 

(iii)  Af  »A, . 

The  case  of  immediate  interest  is  (ii)  in  which  the  free  stream  turbulence 
scale  is  comparable  to  the  boundary  layer  scale.  It  is  felt  that,  while  in  case 
(i)  the  interaction  between  the  free  stream  and  the  boundary  layer  is 
minimal,  case  (iii)  may  require  a  completely  different  approach. 


In  case  (ii),  a  model  can  be  set  up  with  four  layers:  (i)  the  unaffected 
free  stream,  (ii)  the  free  stream  layer  affected  by  interaction,  (iii)  the 
defect  layer  and  (iv)  the  viscous  layer.  Now,  when  asymptotic  matching  of 
the  layers  is  carried  out,  it  is  found  that 

(a)  the  inner  asymptote  of  the  defect  layer  and  the  outer  asymptote  of 
the  viscous  layer  again  become  identical;  and 

(b)  the  magnitude  of  the  Reynolds  stresses  in  the  viscous  layer  becomes 
modified  on  account  of  the  Reynolds  stresses  of  the  free  stream  affecting 
the  defect  layer. 

Result  (a)  is  of  great  significance  because  it  shows  why  a  region  of  log  law 
character  persists  even  with  free  stream  turbulence.  Thus,  in  all  cases,  it 
appears  that  a  viscous  layer  at  the  wall  may  be  matched  with  the  defect 
layer  only  through  the  presence  of  a  log  law  layer.  Furthermore,  as 
Reynolds  number  and  ux  increase,  the  extent  of  the  log  law  region  grows 

with  a  corresponding  (relative)  reduction  of  the  outer  part  of  the  defect 
layer.  In  fact,  the  "defect"  spreads  into  the  free  stream. 

There  are  several  implications  of  the  foregoing,  which  will  be 
discussed. 

1-2.1.  Physical  Structural  Features 

Based  on  the  currently  well-accepted  descriptions  of  the  physical 
structure  of  boundary  layers  in  terms  of  periodic  bursts  and  unstable 
hairpin  vortices,  the  influence  of  free  stream  turbulence  can  be  described 
in  terms  of  the  following. 

(i)  The  perturbed  spanwise  vortex,  which  results  in  a  hairpin  vortex 
with  its  dynamical  features,  also  interacts  with  the  vorticity  of  free  stream 
turbulence.  And, 

(ii)  The  results  is  a  modification  in  the  magnitude  and  distribution  of 
Reynolds  stresses  without  a  modification  of  the  character  of  interactions. 


No  attempt  has  been  made  to-date  on  obtaining  the  modification  to  the 
essential  physical  structure  of  a  boundary  layer  subjected  to  externally 


imposed  vorticity,  particularly  fluctuating  vorticity. 

Writing  the  overall  stress  tensor  in  the  form,  namely 
T ij  =  -  P  8jj  +  2  jx  Sij  =  £  UiUj ,  (2.4) 

one  can  write  the  momentum,  energy  and  vorticity  equations  as  follows. 

UjUij  =  (Tij  /  5) .  j  (2.5) 

P  Uj  (1  Ui  Uj) ,  j  =  (Tij  Uj) .  j  (2  6) 

Uj  (^*) » j  =  -  uj  ®ij  +  Wj  Sij  +  Sij  +  V  ft,jj  (2.7) 


The  mean  vorticity  equation  shows  the  combined  effect  of  fluctuations  of 
vorticity  and  strain.  Whether  additional  mean  or  fluctuating  vorticity  and 
additional  strain  can  give  rise  to  entirely  new  turbulent  processes  is  not 
clear.  However,  currently  available  data  indicate  a  change  in  the  shear 
stress  and  the  outer  intermittency  that  do  not  require  any  new  processes 
other  than  those  currently  accepted  in  a  boundary  layer. 
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1.2.2.  Boundary  Laver  Predictions 

In  the  context  of  LEIH,  it  is  necessary  to  model  eddy-eddy  interactions, 
for  example  in  terms  of  the  skewness  and  the  eddy  viscosity  factors.  In 
predicting  the  (classical)  structure  of  a  boundary  conditions,  the  question 
arises  as  to  the  selection  of  the  two  aforementioned  factors.  The  most 
rational  basis  for  the  selection  of  the  two  factors,  then,  appears  to  be  that 
the  selected  values  of  the  two  factors  must  yield  a  region  of  log  law 
character  adjoining  the  viscous  layer,  both  as  a  necessary  and  a  sufficient 
condition. 

This  has  become  evident  in  predictive  calculations  (reference  30). 
However  the  rationalization  is  based  on  the  results  of  asymptotic  analysis. 

1.3.  STRUCTURAL  SIMILARITY  BETWEEN  VELOCITY  AND 
TEMPERATURE 

The  vast  body  of  literature  on  convective  heat  transfer  appeals  to  the 
possible  analogy  between  heat  and  momentum  transport,  in  terms  of  the 
Prandtl  number,  based  on  the  so-called  Reynolds  analogy,  to  determine 
heat  transfer  to  a  boundary  wall  in  the  fluid.  In  turbulent  flow  this  is 
equivalent  to  drawing  an  analogy  between  eddy  diffusivities  of  heat  and 
momentum.  In  mixing  length  theory,  a  ratio  of  mixing  lengths  for  heat 
and  momentum  is  implied.  The  eddy  diffusivity  for  heat  may  be  defined 

by  —  _ 

Jr  =  -  v0/(dT/0y)  (3.1) 

where  0  is  the  fluctuating  temperature.  The  turbulent  Prandtl  number  is 
then  given  by 


P«T  =  cp  pt/Yr 


(3.2) 


where  Cp  is  the  specific  heat  at  constant  pressure.  The  turbulent  Prandtl 
number  is  assumed  to  be  constant  across  the  boundary  layer,  defined  by 


oT  =  [  uv  /  (8U/3y)]  /  [v0  /  (3T/3y)] 


(3.3) 


where 

U  =  uT  f  (yut/vw  ,  |3) ,  (3.4) 

T  =  Tw  f  (yuT  /  vw  ,  p) ,  (3.5) 

p  being  equal  to  Tc  /  Tw  ;  Tc  is  Squire's  "friction  temperature".  The 
functional  relationships  have  been  specifically  determined  for  various 
classes  of  flows  and  different  parts  of  a  given  turbulent  boundary  layer. 

It  may  be  worth  noting  here  that  there  are  limited,  often  controversial, 
experimental  data  and  largely,  ad  hoc  theories  to  determine  the  influence  of 
molecular  Reynolds  and  Prandtl  numbers  on  Oj  . 

Returning  to  the  consideration  of  different  parts  of  a  boundary  layer, 
there  is  the  outer  boundary  region,  the  interface  that  separates  the 
boundary  layer  and  the  free  stream.  The  free  stream  may  have  additional 
complexities  such  as  pressure  gradient  and  turbulence.  The  interface 
thickness  is  negligibly  small  at  high  Reynolds  numbers,  but  of  the  order  of 
the  sublayer  thickness  or  slightly  larger.  The  pressure  fluctuations  cause 
velocity  fluctuations  over  an  extent  of  the  adjoining  free  stream,  but  there 
are  no  vorticity  fluctuations  in  the  free  stream.  However  the  extent  of  the 
region  over  which  temperature  fluctuations  is  felt  depends  upon  the 
temperature  diffusivity.  The  difference  in  interface  propagation  velocity 
and  velocity  of  entrainment  of  free  stream  fluid  should  also  be  noted.  In 
heat  transfer  problems  one  is  primarily  interested  in  the  entrainment  of 
outer  fluid  into  the  boundary  layer.  The  nature  and  magnitude  of 
entrainment  in  the  presence  of  free  stream  turbulence  is  a  crucial  issue;  the 
influence  of  pressure  fluctuations  must  also  be  included  in  any  calculations. 
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The  other  part  of  a  turbulent  boundary  layer  that  needs  detailed 
consideration  in  heat  transfer  is  the  so-called  inner  layer:  the  combined 
extent  of  the  linear  sublayer  (where  viscous  stresses  are  large  compared  to 
Reynolds  stresses),  the  buffer  layer  (where  viscous  and  Reynolds  stresses 
are  of  the  same  order  of  magnitude)  and  the  logarithmic  law  region  (where 
Reynolds  stresses  are  large  compared  to  viscous  stresses).  As  elsewhere, 
the  logarithmic  law  region  may  be  distinguished  from  the  sublayer,  the 
latter  including  the  buffer  layer  which  must  necessarily  be  thin.  The 
sublayer  is  an  inhomogeneous  region  with  oscillatory  and  time-dependent 
growth  of  streaks  interspersed  between  quiescent  periods.  The  oscillatory 
motion  is  organized  in  several  respects.  Insofar  as  heat  flux  transport  is 
concerned,  it  is  not  clear  if  there  is  a  basic  modification  to  turbulence 
structure  in  the  presence  of  heat  transfer  in  incompressible  flows.  For 
instance,  one  can  pose  the  question  of  whether  there  is  a  difference  in  heat 
transfer  between  the  cases  of  hot  fluid  and  cold  boundary  and  cold  fluid 
and  hot  boundary. 

If  one  considers  the  equation  for  the  rate  of  transport  of  heat,  the 
mechanisms  for  limiting  the  growth  of  the  flux  consist  of  (a)  direct 
dissipation  which  is  small  and  (b)  pressure-temperature  gradient 
correlation.  From  a  structural  point  of  view  the  influence  of  pressure 
fluctuations  must  depend  upon  the  dynamics  of  the  finite  structures  in  the 
boundary  layer,  which,  in  incompressible  flow,  cannot  depend  upon  the 
temperature  of  the  fluid.  Of  course,  the  molecular  viscosity  of  the  fluid  is 
a  function  of  temperature,  but  in  high  Reynolds  number  flows  the  change 
in  viscosity  is  ordinarily  not  effective  in  modifying  diffusion  or  dissipation 
to  any  appreciable  extent. 

On  the  other  hand,  the  entrapment  of  hot  fluid  in  the  spaces  between 
finite  structural  entities  and  its  dynamics  can  influence  heat  transfer 
directly  through  the  influence  of  thermal  conductivity.  It  is  the  influence 
of  such  entrapped  fluid  parcels  (some  cold  and  others  hot)  that  brings  about 
a  difference  between  the  cases  of  relatively  cold  and  hot  boundary  wall. 
Some  aspects  of  entrapment  of  fluid  are  discussed  in  Ref.  31.  However, 
there  is  still  great  uncertainty  in  the  dynamics  of  the  structures;  for 


instance  the  time  scales  associated  with  the  motion  of  the  structures  are  not 
established. 

Both  the  outer  layer  and  the  inner  layer  processes,  of  course,  have  been 
studied  in  the  presence  of  heat  transfer  through  statistical  theory  of 
turbulence  applied  to  space-time  points.  In  this  case,  the  primary  interest  is 
in  the  relation  between  the  spectra  of  velocity  and  temperature. 

In  that  connection,  it  is  of  interest  to  recall  the  following: 

(i)  In  Ref.  32,  the  correlation  UjUjG  has  been  written  as  a  function  of 
q2/e,  ujui,,  uiuj  and  Uj0. 

(ii)  In  Ref.  33,  u$  has  been  modelled  as  a  function^  Q  u,ui 
and  3T/3xi ,  taking  also  into  account  generation  of  1/2  0  and  a  time 
scale  ratio  for  the  decay  of  temperature  and  velocity  fluctuations, 

R  =  [(g2  /  q2)  •  (e  /  e0 )  .  And,  fjnaiiy> 

(iii)  in  Ref.  34,  the  heat  flux  is  modelled  in  the  form 

-v0  =  a0  qV?  (3.6) 

where  ae  is  a  dimensionless  function  of  position  in  the  boundary  layer. 
One  then  needs  a  balance  equation  for  (l/2  02)  which  involves  a  model 
for  the  term  Uj02  .  This  is  written  in  the  following  form. 

-  uj  02  a  02V,  (3.7) 

where  a  so-called  bulk  convection  velocity  is  introduced,  again  as  a 
function  of  position  in  the  boundary  layer.  It  may  be  pointed  out  that  the 
resulting  equations  become  hyperbolic  and  indicate  equilibrium  near  the 
wall.  A  critique  on  this  approximation  may  be  found  in  Ref.  35,  p.  125, 
267  and  298. 
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In  (ii)  of  the  foregoing,  there  arises,  as  pointed  out  earlier,  the  time 
scale  ratio  R,  which  can  also  be  constructed  for  turbulence  intensity  in 
anisotropic  flow.  It  involves  the  variance  of  temperature  fluctuations  and 
the  intensity  of  turbulence.  However,  the  relation  between  Ui0  and  R 
turns  out  to  be  rather  involved.  On  the  other  hand,  considering  a  similar 
relation  for  uluj »  one  can  write 


UjUj  _  q2_ 

vQ  36  Q2  (3.8) 

where  aq  is  a  dimensionless  constant  also.  Similarly,  one  can  consider  two 
bulk  convection  velocities,  one  for  turbulence  intensity  and  another  for 
temperature,  with  two  different  constants  of  proportionality  in  Eq.  3.7.  In 
both  cases,  the  turbulence  intensity  and  temperature  variance  can  be 
expected  to  have  a  relationship  to  the  relevant  spectra. 

1.3.1.  Experimental  Studies  on  Velocity  and  Temperature  Fluctuations 

A  summary  of  experimental  studies  on  velocity  and  temperature 
correlations  and  spectra  in  thermal  boundary  layers  is  presented  in 
Appendix  II  of  this  Report. 

1.3.2.  Modelling  of  Temperature  Fluctuations 

We  concentrate  attention  on  an  incompressible  fluid  flow  of  high 
Reynolds  number,  with  temperature  as  a  passive  contaminant;  the  fluid  is  a 
gas  with  Prandtl  number  approximately  equal  to  unity.  The  mean  heat  flux 
per  unit  area  and  per  unit  time  may  then  be  written  as  follows. 


-7 


dT 


Qj  =  P  cp 


(3.8) 


which  includes  heat  flux  due  to  mean  temperature  gradient  and  that  due  to 
turbulent  transport. 

Introducing  assumptions  on  the  existence  of  a  single  length  scale  and  a 
single  temperature  scale,  one  can  write  the  following: 

0'  =  /3T/ax2  (3.9) 

where  0'  is  the  rms  value  of  0.  Turbulent  heat  transfer  may  then  be 
written  as 


H2  =  -  p  cp  C5  u2  t  3T/0X2,  C5  =  0(1 ) 


(3.10) 


The  justification  for  introducing  a  mixing  length,  t  is  that,  while 
momentum  transfer  k  definitely  not  merely  kinematical  and  depends  upon 
the  dynamics  of  turbulence,  heat  transfer  perhaps  involves  only  a  transport. 


Considering  the  dynamics  of  temperature  fluctuations  for  the  case, 
again,  involving  one  temperature  and  one  length  scale,  it  can  be  argued  that 


0Uj  3T 


dXi 


(3.11) 


In  other  words,  spectral  transfer  of  temperature  fluctuations  towards  the 
small  dissipative  scales  occurs  at  a_rate  governed  by  the  characteristic  time 
of  large  eddies  and  the  amount  of  02  involved. 

It  was  also  shown  by  S.  Corrsin  (Ref.  36)  that  when  the  molecular 
Prandtl  number  is  near  unity,  the  following  can  be  deduced. 

Tle/Tl  (y/v)m  =  1  (3.12) 

under  certain  postulates.  Thus  when  y«  v,  temperature  fluctuations  extend 
all  the  way  up  to  t),  the  Kolmogoroff  scale  for  velocity  fluctuation. 
Molecular  transport  of  heat  is  accelerated  by  straining  motion  in  turbulent 
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flows,  although  mainly  at  small  scales.  Mean  temperature  distribution  is 
strongly  affected  by  turbulence  while  it  is  affected  comparatively  little  by 
molecular  diffusion. 

Finally,  it  is  useful  to  consider  the  spectrum  of  temperature  variance. 
The  three-dimensional  spectrum  can  be  defined  as  follows 

Ee(K)  =  ( f  -1-  4»e  do 

_  JJ  2  (3.13) 

1  ())2  =  f  Eo(K)  dK 

2  Y  J°  (3.14) 

It  represents  the  spectral  density  of  waves  which  have  the  wave  number  of 
magnitude  K  regardless  of  direction  and  is  also  the  spectrum  of 
temperature  variance.  It  may  be  pointed  out  that  <J)e  (K)  is  the  spatial 
spectrum  being  the  Fourier  transform  of  the  spatial  autocorrelation  Re(r). 

Re  (i)  =  exp  (K  ■  r)  exp  (K  •  r)  k  dK 
<t>e  (K)  =  — *— r  If  exp  (-  iK  ■  r)  R*  (r)  dt 

(2ic)3 

As  in  the  case  of  turbulence  kinetic  energy,  one  can  describe  a  cascade 
in  temperature  spectrum.  It  is  noted  that  spectral  flux  of  temperature 
variance  is  the  spectral  equivalent  of  the  following:  u  1  (3T/dX2)  >  where 

u  =  [K2  E(K2)]i/2 

t  =  i/k2  . 


(3.15) 

(3.16) 


The  latter  pertains  to  the  smaller  eddies  that  distort  the  larger  eddies.  The 
spectral  flux  of  temperature  variance  can  then  be  expressed  by 


Te(K)  =  CK2  E0  (KE)l/2 


(3.17) 
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In  the  foregoing,  the  temperature  gradient  associated  with  an  eddy  of  wave 

number  Ki  is  approximated  as  ~  ^  '  It  is  also  useful  to 

point  out  that  the  dissipation  of  temperature  variance  may  be  defined  by 

30  30 
q  =y - 

3Xi  3Xj  (3.18) 

It  can  be  seen  from  the  foregoing  that  the  spectral  features  of 
temperature  fluctuations  depend  upon  the  spectral  features  of  velocity 
fluctuations,  kinematic  viscosity,  diffusivity  of  heat,  Kolmogoroff  scale  in 
the  flow,  dissipation  and  the  inhomogeneities  in  the  mean  strain  in  the  flow. 

In  the  context  of  applying  the  LEIH  to  heat  transfer  problems,  it  can  be 
observed  that  it  is  most  rational  and  also  most  convenient  if  the  spectral 
features  of  velocity  and  temperature  can  be  related  in  anyway.  One  then 
does  not  need  to  introduce  any  assumptions,  for  example,  regarding  the 
variation  of  Prandtl  number  across  a  shear  layer,  which  is  essentially  a 
ratio  of  scales. 

The  question  being  raised  in  the  context  of  the  use  of  the  LEIH  is 
whether  it  is  not  more  satisfactory  to  examine  any  analogy  that  may  exist  in 
the  spectra  of  velocity  and  temperature.  Several  points  of  significance  may 
be  noted  here.  First,  we  referred  earlier  to  the  cascading  of  temperature 
variance;  Eq.  3.17  provides  the  spectral  flux  of  temperature  variance. 
Second,  it  has  been  argued  that  temperature  fluctuations  are  in  a  sense 
transported  by  the  total  velocity.  Temperature  fluctuations  are  related  to 
heat  transfer,  an  "external"  process.  Thus  they  correspond  to  entropy 
fluctuations.  (It  is  possible  to  extend  the  argument  to  the  vorticity 
fluctuations  inherent  to  turbulence,  but  does  not  seem  to  have  been 
pursued).  Third,  the  effect  of  molecular  diffusion  is  felt  only  at  the  small 
scales.  Finally,  while  the  temperature  fluctuation  is  correlated  well  with 
both  ui  and  U2  the  Taylor  spectra  of  ui  and  0  are  different.  In  this 
connection,  one  may  also  note  that  close  to  a  wall,  the  intensities  of  velocity 
fluctuations  are  different,  while  they  are  of  the  same  order  of  magnitude 
elsewhere  in  the  boundary  layer.  The  difference  in  spectra  of  uj  and  0  is 
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in  the  wave  number  range  contributing  most  to  the  variance  of  the  two 
quantities.  At  high  wave  numbers,  an  analogy  can  be  set  up  between  the 
two  spectra.  On  the  other  hand,  at  low  wave  numbers,  where  one  does 
observe  a  difference,  there  is  obvious  interest  in  the  context  of  using  the 
LEIH. 

A  series  of  investigations  undertaken  by  L.  Fulachier,  R.  Dumas,  R. 
Antonia,  Ye.V.  Repik  and  collaborators  on  establishing  relations  between 
the  spectra  of  velocity  and  temperature  (Refs.  37-41)  have  been  described 
in  Appendix  II. 

1.3.3.  Spectra 

In  essence,  it  is  shown  in  the  cited  references  that  a  parameter,  B,  may 
be  constructed  as  follows,  that  is  nearly  constant  across  a  boundary  layer. 

b=£  /  iy 

e2  fy.2  ;)Xl  (3.19) 

The  numerical  value  of  B  from  experiments  is  cited  as  1 .5  in  the  case.  One 
can  compare  this  equation  with  Eq.  3.8. 

The  transfer  terms  have  been  modeled  in  Refs.  X  and  X  by  writing, 

0Uj  =  ae  02  (-  ujuj)] 

It  follows  than  that 

B  =  Prx  ae  /  a1/2  (3.20) 

The  turbulent  Prandtl  number  is  defined  as 
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and  the  argument  is  that  Pry  is  not  necessarily  constant  across  the 
boundary  layer. 

The  authors  cited  above  have  also  obtained  co-spectra  for  u  and  0,  but, 
in  fact,  conclude  that  the  governing  parameter  is  q2  . 

1.3.4.  Correlations  in  the  Near-Wall  Region 

It  may  be  recalled  that  the  physical  case  dealt  with  by  N.  Aubry 
pertains  to  the  structure  of  flow  in  the  near  wall  region  of  a  boundary 
layer.  An  interesting  possibility  in  the  case  of  heat  transfer  is  that  low 
speed  streaks  tend  to  be  heated  when  the  boundary  plate  is  heated  and  high 
speed  streaks  may  consist  of  relatively  cold  fluid  arriving  from  elsewhere. 
Now,  if  the  low  and  the  high  speed  streaks  are  not  dynamically 
interchangeable,  heat  transfer  must  change  in  all  cases,  and  the  structure 
itself  should  change  where  there  are  density  changes. 

In  the  context  of  applying  the  LE1H  to  heat  transfer  problems,  it  is  of 
interest  to  have  some  guidance  from  experiments  on  the  spectra  of  u  and  0. 
Several  observations  in  this  connection  are:  (1)  There  is  near-perfect 
correlation  between  ui  and  0  at  the  wall.  (2)  The  correlation  becomes 
reduced  as  one  moves  out-wards,  and  (3)  at  y+ -  16,  there  is  greater 
correlation  between  spectra  of  the  three  velocity  correlations  than  between 
spectra  of  longitudinal  velocity  fluctuation  and  temperature  fluctuation. 
There  is  close  similarity  between  frequency  spectra  of  velocity  and 
temperature  around  y+  «  2.7. 

There  is  also  some  evidence  from  cross-spectral  data  that  the 
correlation  between  ui  and  0  is  'good'  at  low  wave  numbers.  As  y+ 
increases,  the  wave  number  range  over  which  spectral  coherence  is 
constant  between  ui  and  0  increases;  spectral  coherence  is  defined  by 
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Coh  =  (C§  +  Q2)  <Me  (3.22) 

Co  =  co-spectrum  between  u  and  0 
<{>  =  spectral  density  of  fluctuation  normalized  such  that 

J~  <t>(co)  dco  =  1 

1.3.5.  Spectra  in  the  Logarithmic  Law  Region 

We  may  note  from  Ref.  X  that  it  is  possible  to  conclude  that  the  spectra 
of  u  and  0  are  independent  of  y  in  the  log-law  region.  Further  the  velocity 
spectrum  and  the  temperature  spectrum  are  each  governed  by  the  K*1  law. 

The  main  conclusion  of  this  section  is  that  spectral  similarity  between 

q2  and  02  can  be  invoiced  provided  attention  is  also  paid  to  the  variation 
of  mean  velocity  and  temperature  across  the  boundary  layer. 

Finally,  one  can  introduce  the  concept  of  similarity  in  spectra  and 
determine  mean  temperature  distribution  and  hence  the  skewness  factor  for 
temperature  fluctuation  and  Yr  as  also  the  pdf.  They  can  be  compared 
with  available  experimental  data  on  skewness  factor  and  pdf.  In  the 
alternative,  one  can  formulate  the  temperature  problem  in  the  same  form 
as  the  velocity  problem  and  introduce  a  skewness  factor  for  temperature 
and  Yr  to  obtain  a  solution.  The  resulting  spectra  for  velocity  and 
temperature  can  then  be  compared  with  available  experimental  data.  It 
appears  that  the  latter  procedure  is  likely  to  be  more  difficult  than  the  first. 

1.4.  PROSPECTS  FOR  FURTHER  RESEARCH 

The  main  conclusions  from  the  current  project  are  the  following. 

(i)  As  a  prelude  to  large  eddy  simulation  and  full  Navier-Stokes 
equation  solvers,  the  large  eddy  interaction  hypothesis  provides  a 
relatively  simple  basis  for  establishing  interaction  between  FST  and  BLT  in 
a  fully  developed,  zero-pressure  gradient,  low  speed  TBL.  The  modelling 


involved  is  unambiguous  and  is  related  to  two  factors,  one  in  the  nature  of 
providing  an  excitation  in  the  dynamical  equation  for  the  representative, 
large  eddy  and  the  other  in  the  nature  of  a  damping  or  dissipation  factor. 

(ii)  The  uniqueness  of  a  TBL  rests  on  the  unavoidable  need  for  a 
logarithmic  law-governed  part  of  the  mean  velocity  profile.  The 
emergence  of  such  a  layer  is  often  a  crucial  test  for  the  success  of  a 
prediction  scheme  applied  to  a  TBL.  When  FST  is  present,  it  has  been 
found  numerically  that  a  number  of  solutions  are  feasible  with  given 
boundary  conditions.  The  physically  valid  TBL  solution  is  one  of  them  but 
require  a  major  feature  of  boundary  layer  to  be  incorporated  into  the 
formalism.  Two  possibilities  for  such  a  major  feature  have  been 
considered:  (i)  the  absence  of  a  pressure  gradient  to  first  order  in  the  TBL 
and  (ii)  the  presence  of  a  logarithmic  law-govemed  layer.  Regarding  the 
first,  there  is  still  ambiguity  in  being  able  to  account  for  the  influence  of 
pressure  fluctuations  in  the  free  stream.  On  the  other  hand,  regarding  the 
second,  only  the  existence  of  the  region  is  required  and  not  its  extent, 
which  is  in  fact  a  part  of  the  solution. 

On  the  other  hand,  the  essential  nature  of  interaction  between  1ST  and 
BLT  has  remained  largely  unclear.  Experimental  data  are  still  to  be 
generated  even  in  the  case  of  interaction  between  isotropic  FST  (intensity 
and  scale)  and  fully  developed  TBL  in  respect  of  the  following. 

(i)  Appearance  of  coherent  structures  and  associated  phenomena  in 
streamwise  and  spanwise  directions; 

(ii)  mean  velocity  and  stress  distribution  profiles; 

(iii)  signature  of  pressure  fluctuations  at  the  walls;  and 

(iv)  modifications  in  intermittency,  both  internal  and  at  outer  edge  of 
boundary  layer. 

Such  data  are  urgently  needed  for  elucidation  of  the  effects  of 
interaction.  We  have  emphasized  in  the  past  the  significance  and 


importance  of  conducting  such  investigations  with  free  stream  turbulence 
that  includes  "peaks"  in  the  standard,  uniform  gauze-generated  turbulence. 

In  the  Part  II  of  this  Report,  the  question  of  establishing  free  stream 
turbulence  with  clearly  definable  features  (in  spectrum,  space  and  time)  has 
been  addressed.  It  is  found  that  the  setting  up  of  a  definable  flowfield  is 
not  simple  and  generally  involves  some  ambiguity. 

In  regard  to  modelling  the  interaction,  it  is  useful  to  proceed  as  usual  to 
distinguish  between  the  pointwise  statistical  and  the  large  scale  physical 
structural  approaches.  However,  a  common  link  between  the  two  is  the 
problem  of  determining  what  part  of  FST  is  in  the  nature  of  "inactive" 
motions,  giving  rise  neither  to  stresses  nor  pressure  fluctuations.  The 
problem  is  equivalent  to  that  of  determining  what  part  of  FST  couples  with 
BLT.  In  respect  of  this,  one  may  examine  the  case  of  a  laminar  boundary 
layer  in  the  presence  of  FST.  Experimental  evidence  exists  in  sufficiently 
low  Reynolds  number  flow  wherein  the  mean  velocity  profile,  wall  shear 
stress  and  heat  transfer  are  affected  by  FST  without  the  appearance  of 
transition  to  turbulent  flow.  This  then  represents  the  case  of  a  LBL  that  (a) 
creates  a  balance  between  production  and  dissipation  of  turbulence,  both  of 
small  magnitude  and  (b)  reduces  momentum  deficit  in  itself  through 
acceleration  due  to  turbulent  kinetic  energy  of  the  free  stream.  Detailed 
experimental  data  and  analysis  should  clarify  the  features  of  interaction  in 
this  case,  which  should  also  throw  light  on  the  interaction  between  FST  and 
TBL. 

A  question  of  general  import  pertains  to  the  rational  description  of 
turbulence  in  terms  of  vorticity  fluctuations.  The  large  eddy  interaction 
scheme  has  been  developed  both  for  velocity  fluctuations  and  correlations, 
but  not  for  vorticity  fluctuations.  However,  the  modelling  of  interactions 
between  various  modes  continues  to  be  ad  hoc  in  nature.  In  the  current 
analysis,  for  example,  the  skewness  factor  and  the  turbulent  viscosity  factor 
cannot  be  realized  directly  from  experimental  data. 

A  valuable  direction  for  future  work  on  incorporating  the  physics  of 
turbulence  structure  as  directly  as  possible  into  modelling  appears  to  be  to 


revisit  the  ideas  of  the  single  roller  eddy  due  to  Townsend  and  the 
intermittency  of  turbulent  processes  and  entrainment  due  to  Kovasznay. 
These  ideas  provide  a  framework  for  examining,  especially,  interactive 
fluid  and  thermal  fields.  The  physics  itself  probably  will  remain 
conceptual  in  various  respects. 
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INTRODUCTION 


The  presence  of  free  stream  turbulence  in  practical,  wall-bounded  flows  is  a  common 
occurance.  Whether  the  wall  boundary  layer  is  laminar  or  in  some  state  of  development  us  a 
turbulent  layer,  some  form  of  coupling  between  free  stream  turbulence  (FST)  and  boundary  layer 
rotational  field  can  be  surmized  in  all  cases  except  when  the  scale  and  intensity  of  FST  are  very 
small.  Although  specific  numerical  limits  can  only  apply  in  particular  classes  of  Hows,  one  can 
identify  low,  moderate  and  high  FST  as  shown  schematically  in  figure  1.  In  addition  to  scale  and 
intensity,  FST  may  also  display  inhomogenity,  for  example  a  spectral  gap  or  peakiness  at  various 
wave  numbers,  as  shown  in  figure  2.  The  overall  problem  of  establishing  the  effects  of  FST  must, 
in  general,  include  both  the  states  of  free  stream  and  FST  and  that  of  the  boundary  layer  (laminar, 
traustioning  and  turbulence)  and  also,  of  the  pressure  gradient,  wall  geometry  and  surface 
roughness.  In  the  current  investigation,  we  confine  attention  to  interactions  between  an 
incompressible,  zero  pressure  gradient  turbulent  boundary  layer  and  nearly  isotropic,  moderate 
FST.  Since  an  adiabatic  wall  is  assumed  with  no  heat  transfer,  the  current  case  represents  a  simple 
interaction  between  BLT  and  FST.  Experimental  data  in  such  a  flow  configuration  are  available  in 
Refs.  1  &  2. 

The  problem  of  predicting  the  development  of  a  turbulent  boundary  layer  continues  to  be 
challenging  (Refs.  3-4).  When  FST  is  present,  the  complexity  involved  is  the  interaction  of  a 
random  vortex  fields  with  one  that  has  a  characteristic  structure  in  a  random  field.  The  extent  to 
which  organized  structure  becomes  changed  in  a  TBL  on  account  of  FST  has  not  been  established 
experimentally  to  date.  A  beginning  investigation  has  been  reported  in  Ref.  1  wherein,  through  the 
use  of  a  hot/cold  discimination  technique,  a  burst  frequency  has  been  defined  and  determined  as  the 
average  number  of  hot  intervals  per  unit  time  (when  the  boundary  wall  was  mildly  heated  relative 
to  the  free  stream)  and  the  burst  frequency  has  been  found  to  vary  directly  with  FST  intensity  and 
inversely  with  its  scale.  A  method  of  utilizing  such  data  in  a  predictive  scheme  can  be  visualized, 
but  more  about  this  later.  Meanwhile,  it  is  clear  that  a  method  akin  to  large  eddy  simulation  (l.FS) 
is  probably  required  in  order  to  account  fully  for  the  random  vorlicity  field  interactions. 

In  this  connection,  the  rational  representation  of  turbulence  suggested  by  Lumley  (Refs.  5- 
6)  as  a  means  of  examining  eddy-mean  shear  and  eddy-eddy  interactions,  has  been  chosen  here  for 
application.  The  authors  have  some  experience  with  its  application  to  boundary  layer  flows  (Refs. 
7-10).  I  he  rational  representation  has  been  utilized  to  establish  characteristic  structures  in  a 
channel  flow  (Ref.  1 1),  and  also  to  obtain  the  structure  of  wall  layer  in  a  boundary  layer  (Ref.  12). 
The  Lumley  representation  has  been  utilized  here  entirely  at  the  statistical  level  to  determine  the 
classical  structure  of  TBL  with  FST. 


The  basis  of  statistical  turbulence  representation  is  the  Loeve  orthogonal  decomposition 
theorem  by  means  of  which  the  velocity  fluctuation  is  expressed  in  the  form  ot  generalized 
Fourier  series: 

oo 
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where  n  represents  the  modes,  1,2,3  .....  ^l>  •••  are  random  coelltcients  with  units  of 
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velocity,  uncorrelated  with  one  another.  Yi  are:  structure  functions,  which  are  furthermore 
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assumed  to  be  orthonormal.  The  orthogonality  condition  implies  that  none  of  the  y  i  is 
identically  zero.  It  follows  that 

«n  =  0  (2.1) 

=  ^  ^  Snui  (2.2) 


t>iP)  V,q)  dx  dt  =  5pq 


I  0i  9i  d. 
an  =  j  u  (x, 


t)  <J>in)(x,t)  d  24. dt 


It  is  assumed  that  Ot  j ,  CC2,  OC3  ....  are  ordered  such  that 

A.(1)  >  X(2)  >  l(3)  >  ••••  >  0. 

■  (n) 

Based  on  the  randomness  of  an  and  the  orthogonality  of  9  i  >  one  can  write  for  the 
two-point  correlation, 

00 

Ru  (a,  a' ;  t,  i')  =  X  (tf"’  ♦j"» 


Rij  <t>i")  (x',  t')  dx'  ill'  =  X‘"»  (I.?"  (x  t) 


where  4*  i  are  eigenfunctions  with  as  the  eigenvalues.  For  the  turbulent  kinetic  energy, 
one  can  then  write 
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u2  d&  dt  =  ?i(n) 

^  2  ll=l 


(5) 

(») 

'i 


It  follows  that  X(n)  represents  the  kinetic  energy  content  of  the  entire  flow  associated  with  0  i 
or  the  nth  mode. 

The  uniqueness  of  the  decomposition  is  established  through  calculus  of  variations,  which 
yields  that  each  mode  must  account  for  the  maximum  of  the  energy  available.  One  thus  has  a 
rational,  unique  representation  of  the  velocity  fluctuation  in  terms  of  a  series  of  characteristic 
eddies.  In  a  given  shear  flow,  the  eddies  must  interact  with  one  another  and  with  the  mean  shear 
during  flow  development.  The  immediate  problem  is  partitioning  of  the  energy  among  the  different 
modes.  This  can  be  done  either  by  recourse  to  experimental  data  or  by  truncating  the  nonlinear 
interactions.  On  the  other  hand,  the  large  eddy  interaction  hypothesis  (Llilll)  involves 
consideration  of  only  the  first  mode  over  the  entire  spectrum  and  its  interactions  with  the  mean 
shear  and  all  eddies.  This  is  obviously  a  strong  limitation  of  the  predictive  scheme  presented  here. 

Based  on  Navier-Stokes  equations,  one  can  set  up  a  dynamical  equation  for  the  first  mode 
as  follows. 
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and,  for  incompressible  flow, 
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Ill  Eqns.  (6)  and  (8),  {  }  represents  the  eddy-eddy  interactions.  It  is  essential  to  model  the 

term.  Among  various  possibilities,  the  one  chosen  here  is  as  follows. 


oo  oo 


I  I 
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3xj 
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(9) 


Here  the  first  term  on  the  right  hand  side  involves  a  term  with  a  coefficient  in  the  nature  of  a 
skewness  factor  and  the  second,  one  resembling  eddy  viscosity,  a  type  of  damping  factor.  It  may 
be  stated  at  the  outset  that  the  second  term  with  Vt  is  found  to  be  essential  in  numerical 
predictions.  The  skewness  factor  is  associated  with  the  inhomogeneity  of  the  velocity  field  and  the 
transfer  of  energy  from  regions  of  high  intensity  to  regions  of  lower  intensity. 

F1X)WFIHLC>  MODEL 

The  flowfield  of  interest  is  shown  schematically  in  figure  3  with  four  interactive  regions: 
viscous  layer,  matching  layer,  outer  layer  and  free  stream.  As  stated  earlier,  it  is  assumed  that  the 
no-step  condition  at  the  wall  gives  rise  to  a  boundary  layer  in  the  traditional  sense  for  the  ranges  of 
intensity  and  scale  of  the  FST  present.  Accordingly  the  four  regions  are  described  in  terms  of 
asymptotic  expansions  utilizing  relevant  small  parameters.  There  arise  two  problems. 

First,  the  appearance  of  the  length  scale  of  FST,  Le,  in  addition  to  the  inviscid  length  scale 
of  the  flow,  1:  there  is  no  difficulty  in  constructing  the  required  expansion  parameters  utilizing  Lc, 
free  stream  velocity,  Uc,  boundary  layer  thickness  6  and  the  outer  region  length  scale,  A  . 
However,  the  describing  equations  involve  no  physical  influence  of  the  FST  length  scale. 
Therefore,  we  have  proceeded  as  follows:  a  value  of  60.995  *s  chosen  and  the  dissipation  length 
scale,  Le  ,  is  determined  corresponding  to  the  given  value  of  FST  intensity;  and  a  boundary 
condition  is  introduced  at  the  edge  of  the  boundary  layer  requiring  the  dissipation  length  scales  in 
the  boundary  layer  and  the  free  stream  to  be  equal  at  the  interface.  In  other  words,  that  solution  is 
found  which  yields  a  dissipation  length  scale  within  the  boundary  layer  equal  to  that  in  the  free 
stream.  The  matched  asymptotic  expansion-solution  procedure  is  otherwise  based  on  an  extension 
of  the  method  of  Refs.  13-14.  including  the  presence  of  FST.  The  other  problem  is  concerned 
with  the  matching  layer.  In  Refs.  1-2  and  15-16  there  is  considerable  evidence  for  the  existence 
and  persistence  of  the  log  law  character  of  the  region  and  also  for  the  transverse  growth  of  the 
region  with  FST  intensity.  In  Ref.  4  for  the  case  of  an  ordinary  TBL,  one  without  FST,  an  attempt 


was  successful  in  predicting  the  existence  of  a  constant  stress  layer  corresponding  to  the  log  law 
region.  A  coordinate  transformation  and  some  related  assumptions  are  difficult  to  understand.  In 
Ref.  17  the  log  law  character  of  mean  flow  variation  was  prescribed  in  dealing  witli  an  ordinary 
TBL.  When  FST  is  present,  as  stated  above,  there  is  a  need  to  match  BLT  with  FST  quantities.  In 
order  to  assure  a  turbulent  boundary  layer  type  solution,  it  is  then  found  to  be  necessary  to  invoke 
the  boundary  layer  character  in  some  way.  In  the  current  formulation  it  is  assumed  that  the  mean 
flow  in  the  matching  layer  is  governed  by  the  log  law.  Prescribing  this  does  not,  of  course, 
specify  the  extent  of  the  log  law  region  or  that  of  the  defect  region. 

TEST  FLOW  CONFIGURATION 

The  flow  configuration  chosen  is  described  in  detail  in  Refs.  1  and  2.  A  wind  tunnel  with  a 
speed  of  about  16.5  m/sec.  was  utilized  to  examine  a  flat  plate  boundary  layer  while  the  FST 
intensity  and  length  scales  were  varied  by  changing  the  upstream  grids.  The  maximum  intensity  at 
any  location  has  been  reported  to  be  about  6  per  cent,  all  such  locations  yielding  Ue  0  /  v  greater 
than  2,000.  The  investigation  covered  various  aspects,  but  the  details  regarding  the  few  test  cases 
chosen  for  prediction  are  summarized  in  figures  4-6  and  Table  I.  It  is  of  particular  interest  to 
observe  (a)  the  simultaneous  influence  of  intensity  and  length  scale  of  FST,  (b)  the  persistence  and 
enlargement  of  the  log  law  character  of  the  matching  layer  in  the  boundary  layer  region,  and  (c)  the 
magnitude  of  Le  /  80.995  as  a  function  of  Lc  ,  intensity  and  Reynolds  number. 

FLOWFIELD  PREDICTION  PROCEDURE 

The  basic  numerical  method  is  a  fractional  method  (Refs.  19,  20)  applied  to  equations  for 
the  mean  flow  and  also  for  the  large  eddy  functions  <j>j  .  The  describing  equations  for 
incompressible,  time-dependent  flow  are  written  in  conservative  form.  The  pressure  field  is 
determined  from  the  Poisson  equation.  The  crossflow  direction  z  is  transformed  into  wave 
number  space  k.  The  convection  terms  are  transformed  with  second  order  Adams-Bashforth 
differencing,  and  the  viscous  terms  with  implicit  second  order  Crank-Nicolson  method. 
Regarding  the  grid,  a  staggered  grid  is  employed  for  discretization:  the  u  -  component  variables 
are  placed  on  the  east-west  side  of  the  grid  cell,  the  v  -  component  ones  on  the  south-north  side  and 
the  w  -  component  variables  and  pressure  are  located  at  the  cell  centers.  The  resulting  Poisson 
equations  are  solved  implicitly  with  an  elliptic  procedure  without  resorting  to  an  approximate 
factorization  technique. 

The  mean  flow  equations  are  decoupled  from  the  large  eddy  equations  during  integration  in 
time  by  allowing  the  turbulent  transport  terms  in  the  mean  flow  equations  to  lag  by  the  time  step 
At  when  the  mean  velocities  are  being  solved  at  t=  (tu  +  At ).  Meanwhile  the  large  eddy  field, 
interacting  with  the  current  mean  flow  field,  is  implicitly  updated  for  t  =  (tu  +  At  ). 


Other  aspects  of  the  numerical  scheme  are  described  in  Appendix  I.  In  general, 
approximately  3(W)  iterations  are  required  in  time  at  each  point  in  the  flowfield  for  the  chosen  values 
of  S  and  vt .  Then  three  sets  of  matching  conditions  have  to  be  introduced  along  the  transverse 
direction  invoking  the  log  law  character  of  mean  flow  in  the  matching  layer.  Finally,  the  calculated 
flowfield  has  to  be  iterated  with  respect  to  matching  boundaries  till  the  boundary  conditions  are 
satisfied.  The  solution  so  obtained  applies  to  the  free  stream  turbulence  case  of  given  intensity  and 
a  dissipation  length  scale  as  obtained  at  the  outer  edge  of  the  boundary  layer. 

Preliminary  Investigation 

In  order  to  test  the  procedure,  the  standard  equilibrium  boundary  layer  case,  the  one 
equivalent  to  the  Falkner-Skan  family  of  laminary  boundary  layers,  (Ref.  18),  was  computed  for 
Ree  of  2,500.  The  convergence  of  the  solution  as  well  as  the  method  of  choosing  the  values  of  S 
and  vt  were  found  satisfactory.  It  may  be  noted  that  since  there  is  no  free  stream  turbulence  in  this 
case,  it  is  adequate  to  demand  convergence  of  the  solution. 

PREDICTIONS  AND  DISCUSSION 

A  comparison  of  experimental  data  and  predictions  is  presented  in  respect  of  the  following 
in  various  test  cases:  mean  flow  in  figures  7-11,  shear  stress  variations  in  figures  12-14  and 
intensity  variations  in  figures  15-18.  Some  difficulties  in  predictions  that  arose  in  the  beginning  at 
the  matching  stations  were  corrected  subsequently  by  an  adjustment  of  the  computational  grid.  It  is 
worth  mentioning  that  it  was  necessary  to  utilize  two  different  values  for  the  skewness  factor  in 
the  x  and  y  directions  in  any  given  case;  however,  the  values  were  applicable  throughout  the 
boundary  layer.  On  the  other  hand,  the  values  of  Vj  were  generally  immaterial  so  long  as  a  non¬ 
zero  value  was  included. 

Discussion 

The  investigation  pertains  to  a  turbulent  boundary  in  a  free  stream  of  moderate  turbulent 
intensity  and  scale.  The  interaction  region  extends  over  the  scale  of  FST. 

The  predictive  scheme  rests  on  a  number  of  premises.  First,  in  order  to  assure  a  turbulent 
boundary  layer  in  the  presence  of  FST,  the  matching  layer  between  the  viscous  sublayer  and  the 
outer  layer  needs  to  be  specified  as  governed  by  the  Log  Law.  This  is  based  on  the  fact  that  other 
types  of  boundary  layers  are  also  feasible  in  a  FST  environment,  and  also  on  the  lack  of  a 
mechanistic  model  of  interaction  between  FST  and  BLT.  Second,  the  outer  edge  of  the  boundary 
layer  is  treated  as  a  parameter  for  establishing  matching  of  turbulence  dissipation  scales.  No 
distinction  is  made  between  the  locations  where  the  values  of  turbulence  intensities  and  mean  flow 
are  matched,  although  the  presence  of  the  wall  must  affect  the  y  -  component  of  intensity.  In  the 
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actual  computational  work,  experimental  data  were  utilized  extensively  in  reducing  the  irinl-and- 
error  procedure  required  in  matching  free  stream  and  boundary  layer  Hows.  Third,  the  application 
of  LEIH  is  based  on  the  representation  of  turbulence  by  the  first  mode  following  orthonomal 
orthogonal  decomposition.  Furthermore,  there  is  no  physical  significance  to  the  definition  of 
skewness  and  eddy  viscosity  factors. 

It  is  gratifying  that  a  choice  of  skewness  and  eddy  viscosity  factors  is  adequate  for 
reproducing  the  main  turbulence  quantities  as  measured  in  an  experiment.  In  the  predictions,  it 
has  (obviously)  been  possible  to  obtain  spectra;  details  are  provided  in  Ref.  21.  The  main 
observations  from  the  spectra  are  (i)  the  contributions  from  small  wave  numbers  becomes 
increased  when  FST  is  present,  more  than  those  from  large  wave  numbers;  and  (ii)  the  changes  for 
v2  are  smaller  than  the  changes  for  u2. 
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Table  I. 


VALUES  OF  SF  AND  vtf  UTILIZED 


Spj 

sFj 

Vjpj 

V  1 1  j 

Case  1 

0.8 

0.82 

0.1  1 

0.10 

Case  2 

0.92 

0.92 

0.1  1 

0.10 

Case  3 

0.91 

0.93 

0.1 1 

0.10 

Case  4 

0.90 

0.93 

0.1 1 

0.10 

Case  5 

0.93 

0.91 

0.1 1 

0.10 
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Appendix  I 


The  mean  flow  equations  are  written  a  conservative  form  as  follows: 


dUj 

at 


a2  u, 

ax? 


(A.l) 


with 


UjUj  =  <t>i  4>j 


(A. 2) 


The  large-eddy  functions  <))j  are  governed  by  the  dynamical  equation  in  physical  (x,y,z)  space  as 
follows. 


a<t>i 

at 


+  -  (^i  Uj)  +  — —  (<))i  Uj)  + 

dxj  Dxj 


S 


dn 

= - +  (V  +  vt) 

3x; 


ax? 


(A. 3) 


In  the  following,  the  solution  procedure  is  illustrated  only  for  the  x-component  of  <h  .  The 
same  procedure  can  be  extended  to  other  components  of  <t>j . 

The  flow  is  assumed  to  be  homogeneous  in  the  spanwise  z  -  coordinate  direction.  One  can 
then  introduce  Fourier  expansions  for  velocity  and  pressure  as  follows. 


<J>i  (x.y.z.t)  =  £  <t>i  (x,y,k,t)  eikz 

k  =  -  oo 


7t  (x,y,z,t)  =  X  71  (x.y.k.O  eikz 

k  =  -  oo 


(A  .4) 


(A. 5) 


where  i=  V-  1  .  Applying  these  to  Eq.  (A. s),  one^obtains  a  complex  equation  for  4>i  (x,y,k,t) 
as  shown  below. 
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—  +  2  —  (u  4>i)  +—  (v  $,)  +  —  (u  (J>2 

dt  dx  dy  dy 


-  (  d  1 

+  ikU<j>3  +  S  •  F  —  (<>i <t>j ) 

l3xj 


—  +  (v  +  vj  (—  +  —  -  k2] 
dx  \dx2  dy2 


(A. 6) 


where  the  F  -  term  is  equivalent  to 


OO  / 

X  -k  (k')  01  (k  -  k')]  +  —  [<t>i  (k')  <J>2  (k-k') 
k'  =  -  oo  l 3x  Dy 


+  i  k  <)>,  (k')  4>3  (k  -  k')J 

Thus  infomiation  between  a  large  and  a  small  wave  number  is  exchanged  through  the  convolution 

/\ 

term.  The  solution  for  <j>i  can  then  be  obtained  at  time  t  =  to  +  At  utilizing  the  values  of  mean 
velocities  U  and  V  which  have  been  determined  corresponding  to  that  instant  but,  as  stated  under 
Numerical  Procedure,  with  the  Reynolds  stresses  for  t  =  t0. 

A 

* 

The  first  step  in  solving  Fq.  (A. 6)  is  to  determine  an  intermediate  function  0*  from  the 
following  equation. 


01-01  .  ,,  1^01  ctyi  dU  ^  2  \ 

-  +  U  -  +  V  -  +  -  (pi  +  -  <{>2 

At  dx  dy  dx  dx  / 


+  1  {3  S  F  4,  <J)j)n  -  S  F  | 


-  I  1  +  Vi  iL.  +  iL  _  k2 

^  tdx2  dy2 


(m 


(A. 7) 
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It  is  subject  to  the  conditions 

n 

(A. 8) 


*  "'ii+i 

01  =  <t>i 


+  (At) 


^r1 


\d*i 


^+1 
“  0i 


+  (At) 


c)  7t 

{d*t 


vn+ 1 


Equation  (A. 8)  assumes  that  k  ~  n  at  the  boundary  and  is  a  second-order 

A 

* 

approximation  in  time,  O  (At2)  (Ref.  20).  Introduction  of  in  the  viscous  term  in  Eq.  (A. 8) 

A.  <V, 

stabilizes  the  numerical  solution  for  (t>1  .  Once  ^i  is  solved  from  the  elliptic  equation  (A. 8), 
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then  the  pressure  spectra  ji  is  obtained  from 
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The  intensities  and  Reynolds  stress  components  Ujiij  are  then  calculated  from  (j>j  , 
integrating  over  a  chosen  wave  number  range  as  follows. 
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where  *  on  denotes  the  complex  conjugate.  It  is  found  that  a  set  of  discrete  wave  numbers 
encompassing  the  desired  range  is  adequate  in  practice. 
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Figure  4.  Experimental  set-up. 
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Figure  15.  Direct  stress  prediction:  case  I 
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An  attempt  is  made  to  illustrate  a  relation  between  large  eddy  interaction 
hypothesis  and  rapid  distortion  in  the  context  of  a  shockwave  and  an  expansion 
interacting  with  turbulence  as  depicted  in  the  experimental  studies  of  J.  F.  Debieve 
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ABSTRACT 


An  attempt  is  made  to  illustrate  a  relation  between  large  eddy 
interaction  hypothesis  and  rapid  distortion  approximation  in  the  context 
of  expansion  of  turbulent  supersonic  flow,  J.F.  Debieve  et  al.  1982  based 
on  J.P.  Dussange,  1981.  Three  cases  are  compared  in  relation  to 
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wall.  The  differences  between  the  three  cases  are  illustrated  in  terms  of 
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pressure-strain  correlations.  Similar  analysis  is  underway  for  velocity- 
temperature  correlations. 
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•  BACKGROUND  (J.P.  DUSSAUGE,  1981). 

•  FORMULATION  IN  MASS-WEIGHTED  VARIABLES 

•  INTRODUCTION  OF  VARIABLE  THAT  ACCOUNTS  FOR 

MEAN  DILATION. 

•  RAPID  DISTORTION 

•  DENSITY  VARIATIONS  FROM  MEAN  MOTION 

•  SOLENOID AL  VELOCITY  FLUCTUATIONS 

•  SIMILARITY  BETWEEN  U,  AND  0 

•  PRESSURE  STRAIN  WITH 

•  INCOMPRESSIBLE  p'  AND 

•  Djj  INSTEAD  OF  dn-Jdx j 

•  MAIN  RESULTS 

•  TURBULENCE  INCREASES  ACROSS  SHOCKWAVE 

•  AT  HIGH  WAVE  NUMBER 

•  TURBULENCE  DECREASES  ACROSS  EXPANSION 

•  WITH  UNEXPLAINED  BEHAVIOR  AT  HIGH  WAVE 
NUMBERS 

•  INTEGRAL  SCALE  DEPENDS  STRONGLY  ON 
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•  LARGE  EDDY  INTERACTION  HYPOTHESIS 
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under  the  assumption  of  linearization  and  neglect  of  ff. 
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•  CALCULATION  PROCEDURE 


•  REGION  OF  INTEREST 

•  SUPERSONIC  FLOW 

•  BOUNDARY  LAYER  WAKE 

•  OUTER  TURBULENT  FLOW 

•  FIND  (f  FROM  p'  and  0 
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•  ESTABLISH  SPECTRAL  CONTENT 

•  COMPARE 

•  RAPID  DISTORTION  OVER  CURVED  WALL 

•  NORMAL  DISTORTION  OVER  CURVED  WALL 

•  INCOMPRESSIBLE  FLOW  OVER  CURVED  WALL 

•  TWO  SETS  OF  RESULTS 

•  LE1H  UNDER  RAPID  DISTORTION  APPROXIMATION 
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•  MAIN  RESULTS 

•  PRESSURE-STRAIN  CORRELATIONS  SEEM  TO 
PRESENT  SOME  AMBIGUITIES  IN  RAPID  DISTORTION 
APPROXIMATION  COMPARED  TO  TURBULENCE 
INTENSITY. 

•  PRESSURE-STRAIN  CORRELATION  ESTIMATES 
SHOW  THE  EFFECTS  OF  INACTIVE  MOTION 
WHICH  BECOME  DECREASED  DURING  EXPANSION. 
HOWEVER,  DISSIPATION  SEEMS  TO  INCREASE. 

•  RAPID  DISTORTION  APPROXIMATION  YIELDS 
RESULTS  OF  DOUBTFUL  VALIDITY  EXCEPT 
AT  VERY  LOW  WAVE  NUMBERS. 

•  THE  LEIH  IS  SUITABLE  FOR  EVALUATING 
COMPRESSIBILITY  EFFECTS  IN  THE  PRESENCE 
OF  FREE  STREAM  TURBULENCE. 

•  THE  LERI  MODEL  VERIFIES  EXPERIMENTALLY 
OBSERVED  DIFFERENCES  BETWEEN  INNER, 

OUTER  AND  FREESTREAM  REGIONS. 

•  CORNER  EXPANSION  YIELDS  RESULTS 
SIMILAR  TO  RELAXING  FLOWS. 

•  CURVATURE  DOMINATES  OVER  EXPANSION 
EFFECTS. 
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II. 1 


Summary  of  Imperial  College  work  on  free-stream  turbulence 
(ref.  no.  989AQ/AB)  under  contract  with  Purdue  University. 

P.  Bradshaw,  Principal  Investigator 
Present  address:  Mechanical  Engineering  Dept,  Stanford  University 

1  August  1988  (updated  29  August  1990). 


1 •  Introduction 

The  main  object  of  this  work  was  to  study,  experimentally, 
the  effects  of  anisotropic  free-stream  turbulence  on  heat 
transfer  in  low-speed  turbulent  boundary  layers,  with 
particular  reference  to  turbomachinery  applications.  Typical 
free-stream  turbulence  intensities  were  of  the  order  of  5  percent  of 
the  free-stream  velocity  -  that  is,  typical  of  turbomachine  blade 
wakes  rather  than  the  efflux  from  combustion  chambers. 

The  experiments  were  done  in  three  stages:  first,  basic  measurements 
with  isotropic  free-stream  turbulence,  to  compare  with  existing 
data;  second,  studies  of  the  effect  of  turbulence  generated  by 
oscillating  spanwise  rods  (a  traditional  method  of  simulating  turbo¬ 
machine  turbulence);  and  third  -  following  the  abandonment  of 
Stage  Two,  for  reasons  to  be  explained  below  -  measurements  of  the 
effect  of  the  wake  of  a  stationary  rod  slightly  inclined  to  the 
spanwise  direction,  to  simulate  non-radial  blades.  The  data  were 
intended  to  feed  directly  into  the  large-eddy  interaction  model  of 
turbulence  being  developed  at  Purdue  University.  The  work  at 
Imperial  College  terminated  when  the  Principal  Investigator  left  in 
August  1988:  fortunately  this  coincided  with  a  natural  break  in  the 
program,  and  the  work  continued  at  Stanford  University  in  a  slightly 
different  direction.  The  main  publications  to  date  have  been  by 
Baskaran  &  Bradshaw  (1989)  and  Baskaran,  Abdellatif  &  Bradshaw  (1989). 


2.  Equipment. and  Techniques 

The  measurements  were  all  made  in  the  Imperial  College  Aeronautics 
Department  30  in  x  30  in  closed-circuit  wind  tunnel,  whose  natural 
turbulence  level  is  approximately  0.03  percent,  two  orders  of 
magnitude  below  the  level  artificially  generated  in  the  present  study. 
The  test  boundary  layer  was  genrated  on  a  flat  plate  8  ft  in  chord, 
with  the  ogive-shaped  leading  edge  developed  by  Hancock  to  avoid 
leading-edge  separation  at  rms  v-component  turbulence  intensities 
up  to  approx.  5  percent.  A  trailing-edge  flap  was  fitted  to  adjust 
the  circulation  around  the  plate  to  ensure  symmetrical  flow  over 
the  leading  edge.  The  plate  was  generally  similar  to  that  used 
by  Simonich  &  Bradshaw:  strips  of  steel  shim  stock  1  in  wide  and 
approx.  0.001  in  thick  ran  fore  and  aft,  sandwiched  between  two 
sheets  of  plywood  approx.  1/4  in  thick  and  electrically  connected 
in  series.  When  the  plate  was  heated  (with  an  a.c.  voltage  of  about 
100V,  giving  a  current  of  13A),  half  the  heat  was  transferred  through 
the  top  surface  and  half  through  the  bottom  surface,  so  that  the 


(spatially-uniform)  rate  of  heat  transfer  could  be  deduced  directly 
from  the  electrical  power  input,  with  no  correction  for  losses  except 
those  due  to  radiation.  The  surface  temperature  was  measured  with 
r.-vi rc  thermocouples,  attached  to  the  surface  with  varnish  which 
reduced  the  effective  roughness  well  below  the  critical.  The  cold 
junctions  of  the  thermocouples  were  mounted  on  top  of  plexiglass 
posts,  projecting  from  the  plate  surface  into  the  free  stream,  well 
to  one  side  of  the  centre  line  on  which  the  measurements  were  made. 
Surface  shear  stress  was  deduced  either  from  Preston  tube  measure¬ 
ments  or  by  fitting  mean-velocity  profiles  to  the  law  of  the  wall. 
Standard  constant-temperature  hot-wire  equipment  was  used.  More 
details  are  given  by  Baskaran  &  Bradshaw  (1989). 


3 .  Results 

3.1  First  Stage  -  isotropic  free-stream  turbulence 

Measurements  of  heat  transfer  and  skin  friction 
on  flat  plates  behind  conventional  square-mesh  grids,  as 
used  at  Imperial  College  by  Simonich  (Simonich  &  Bradshaw  1978) 
and  by  Hancock  (Hancock  &  Bradshaw  1989  and  papers  cited  therein), 
were  made  as  checks  of  techniques  and  of  Simonich 's  early  data, 
obtained  before  Hancock's  demonstration  of  large  effects  of  free- 
stream  turbulence  length  scale  on  the  skin  friction  for  a  given 
free-stream  turbulence  intensity.  The  results  generally  confirm 
Simonich's  but  show  that  Hancock's  correlation  for  the  combined 
effects  of  free-stream  intensity  and  length  scale  does  not  work 
very  well  for  heat  transfer  (i.e.  Stanton  number ) although  it  is  an 
excellent  collapse  of  many  workers'  skin-friction  results.  The 
discrepancy  is  not  merely  a  departure  from  Hancock's  correlation 
curve:  judging  by  the  present  data,  Hancock's  combined  length-scale 
and  intensity  parameter  does  not  provide  a  unique  collapse  of 
fractional  changes  in  Stanton  number  as  it  does  for  fractional 
changes  in  skin-friction  coefficient.  The  variations  in  Reynolds 
analogy  factor  implied  by  the  present  measurements  are  larger  than 
in  the  work  of  Blair  (1983). 

At  the  end  of  the  program  a  complete  set  of  u-component  spectra 
was  measured,  both  throughout  the  boundary  layer  and  in  the  free 
stream,  to  define  length  scale  behavior  in  more  detail  (the  length 
scale  used  for  most  of  the  data  correlations  being,  for  simplicity, 
the  "dissipation  length  parameter"  derived  from  estimates  of 
dissipation  deduced  from  free-stream  turbulence  decay  rate). 

3.2  Second  stage  -  wakes  of  spanwise  rods 

Measurements  were  made  of  the  skin  friction  on  a  flat  plate  behind 
spanwise  circular  cylinders,  either  stationary  or  oscillating  in  the 
direction  normal  to  the  plate.  Arrays  of  cylinders  in  the  form  of  a 
spoked  wheel,  rotating  about  an  axis  in  the  stream  direction,  have 
been  used  by  many  workers  to  simulate  the  effect  of  rotor  wakes  on 
a  stator,  or  vice  versa.  A  mechanism  to  oscillate  horizontal  rods 
in  the  vertical  direction  was  already  available  in  the  30  in  tunnel. 
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for  use  in  studies  of  vortex  streets  behind  oscillating  bluff 
bodies.  We  therefore  elected  to  use  a  single  rod,  oscillating 
through  a  sufficient  amplitude  in  the  vertical  direction  that  the 
vertical  velocity  of  translation  of  the  wake  during  its  interaction 
with  the  plate  boundary  layer  was  essentially  constant. 

Initial  tests  were  made  with  the  rod  stationary,  with  its  axis 
in  the  plane  of  the  plate  so  that  the  wake  met  the  plate  leading 
edge  symmetrically.  The  results  showed  reduced  skin  friction, 
which  was  unexpected  in  view  of  the  increases  reported  by  users  of 
spoked-wheel  rod  arrays.  The  explanation  was  qualitatively  simple, 
if  quantitatively  surprising:  the  effective  free-stream  velocity 
over  the  plate  is  reduced  because  of  the  reduced  velocity  in  the 
wake,  and  evidently  this  outweighs  the  effect  of  the  turbulence  in 
the  wake.  This  result  is  not  special  to  the  present  generator  and 
would  probably  be  found  with  any  stationary  wake  generator.  Reduced 
skin  friction  was  also  found  when  the  rod  was  oscillated  vertically. 
However  a  complete  row  of  rods  or  blades,  rotating  or  not,  would  tend 
to  increase  skin  friction.  This  is  because  skin  friction  depends, 
roughly,  on  (free  stream  velocity)2,  and  the  increase  of  axial 
speed  between  the  wakes,  required  to  satisfy  continuity  (constant 
volume  flow  rate  in  the  present  low-speed  experiment)  would  outweigh 
the  effect  of  the  velocity  defect  in  the  blade  wakes. 

These  preliminary  results  rather  spectacularly  exposed  the  lack  of 
relevance  of  the  many  previous  studies  of  moving  rod  wakes,  and  we 
decided  to  abandon  this  kind  of  turbulence  generator.  The 
continuation  of  this  work  at  Stanford  University,  under  AFOSR 
grants  AFOSR-86-0073  and  89-0246,  uses  a  two-part  grid  to 
produce  a  turbulent  wake  without  a  mean  velocity  gradient,  so 
that  the  effect  of  turbulence  as  such  can  be  studied  in  the  absence 
of  mean-velocity  variations. 

3.3  Third  stage  -  wake  of  a  stationary  inclined  rod 

An  important  kind  of  wake-blade  interaction  occurs  when  the  wake 
is  not  parallel  to  the  blade  leading  edge.  An  upstream  rod 
IS  a  useful  simulation  of  this  kind  of  wake,  if  "rolled"  (in  the 
aeronautical  sense  of  the  term)  at  an  angle  to  the  span  of  the  flat 
plate  representing  the  downstream  blade.  Almost  certainly,  spanwise 
anisotropy  /  inhomogeneity  is  more  important  than  variations  normal 
to  the  plane  of  the  blade,  so  this  is  a  key  configuration  to  study. 

Up  to  the  time  of  closure  of  the  Imperial  College  work,  measurements 
of  all  mean-velocity  components,  all  Reynolds  stresses  and  all 
triple  products  of  velocity  fluctuations  had  been  made  for  rod  roll 
angles  of  0,  10  and  20  deg.,  and  for  a  flat-plate  boundary  layer 
in  the  absence  of  a  rod  (as  control).  A  full  map  of  the 
Stanton  number  variation  over  the  plate  surface  was  also 
obtained  for  all  three  rod  angles.  An  obvious  feature 
is  that  although  spanwise  variations  are  quite  sharp,  too  sharp 
for  the  3D  "boundary  layer"  approximation  to  be  applied,  the 
W  (spanwise)  component  of  mean  velocity  and  the  uw  and  vw 
Reynolds  stresses  were  fairly  small,  suggesting  that  quasi-2D 
"strip  theory"  calculations  might  be  adequate  in  practice.  We 
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have  not  yet  used  these  data  for  test  cases  for  calculation 
methods  but  hope  that  this  could  be  done  in  future. 


4.  Conclusions 

Although  the  program  at  Imperial  College  was  cut  short,  it  yielded 
several  useful  pieces  of  information  as  well  as  providing  data  and 
insight  for  the  modelling  program.  The  simple  configuration  chosen  to 
represent  a  turbomachine  blade  surface  (a  flat  plate  in  low-speed 
flow)  removed  subsidiary  phenomena  like  compressibility  and  pressure 
gradient,  allowing  clearer  conclusions  to  be  drawn,  (i)  The  failure 
of  the  heat-transfer  results  for  isotropic  free-stream  turbulence  to 
collapse  on  the  combined  length-scale/intensity  variable  chosen  by 
Hancock  to  collapse  skin-friction  data  suggests  that  there  are 
significant  physical  differences  in  the  effects  of  free-stream 
turbulence  on  heat  transfer  and  on  skin  friction  (momentum  transfer) . 
(ii)  The  decrease  in  skin  friction  found  on  a  flat  plate  in  the  wake 
of  a  single,  stationary  simulated  blade  exposed  the  misleading  nature 
of  previous  workers'  demonstrations  of  increased  skin  friction  behind 
arrays  of  moving  "blades":  there  is  a  simple  explanation  based  on 
continuity  arguments,  (iii)  Measurements  on  a  flat  plate  downstream 
of  a  rod  inclined  in  front  view,  simulating  a  non-radial  blade, 
showed  that  although  spanwise  gradients  were  quite  large,  the  spanwise 
mean  velocity  and  the  shear  stresses  with  spanwise  components  were 
small:  this  may  legitimise  the  use  of  2-D  "strip  theory"  methods  for 
predicting  boundary  layers  on  such  blades. 
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An  experimental  study  into  the  wake-boundary  layer  interaction 
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Synopsis 

Mean  flow  and  turbulence  measurements  were  conducted  in  a  flat  plate  turbu¬ 
lent  boundary  layer  merging  with  the  wake  of  a  circular  cylinder  inclined  in  the 
cross-stream  plane,  simulating,  to  a  first  approximation,  the  wake  of  an  axial  tur¬ 
bomachine  turbine  blade  interacting  with  the  boundary  layer  over  another  blade 
downstream.  Experiments  were  conducted  for  angles  of  inclination  of  the  wake, 
with  respect  to  the  plane  of  the  flat  palte,  of  0°,10°  and  20°.  The  results  demon¬ 
strate  that  the  primary  effect  of  the  interaction  is  simply  to  reduce  the  skin  friction 
due  to  the  reduced  velocity  in  the  wake  of  the  cylinder.  The  wake  is  bodily  shifted 
towards  the  local  free  stream  by  the  growing  boundary  layer.  The  turbulence  of 
the  wake,  regarded  as  “free  stream  turbulence”,  appears  to  have  little  effect  on  the 
skin  friction. 


t  Thermosci.  Divn.,  ME  Dept.,  Stanford  University,  Stanford,  CA  04305. 
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Introduction 

This  paper  deals  is  a  part  of  the  investigation  into  “complex  ”  flows  at  Imperial 
College,  it  is  one  of  several  papers  on  the  interaction  of  a  simple  shear  layer  with 
another  turbulent  field,  in  this  case  a  flat  plate  boundary  layer  interacting  with 
the  wake  generated  by  a  circular  cylinder.  The  present  work  differs  from  earlier 
studies  on  interaction  between  circular  cylinder  wakes  and  boundary  layers  (see  for 
example  Marumo,  Suzuki  ic  Sato  1978)  in  that  the  circular  cylinder  was  positioned 
well  upstream  of  the  leading  edge  of  the  flat  plate  on  which  the  boundary  layer 
grows,  whereas  in  the  earlier  investigations  the  cylinder  was  held  inside  the  turbulent 
boundary  layer  and  parallel  to  the  plate.  The  present  flow  configuration  in  effect 
idealises  the  flow  in  a  turbomachine,  where  the  wake  of  a  blade  interacts  with  the 
boundary  layer  over  another  blade  downstream. 

Since  the  leading  edge  and  the  trailing  edge  lines  of  a  blade  are  in  general  not 
parallel,  due  to  blade  twist,  the  wake  meets  the  span  of  the  downstream  blade  at  an 
angle.  As  a  consequence,  the  skin  friction  drag  and  the  surface  heat  transfer  rates 
over  the  blade,  and  hence  the  overall  performance  of  the  turbomachine,  is  expected 
to  be  affected,  in  a  more  complicated  way  than  if  the  blades  were  parallel.  There 
appear  to  be  practically  no  published  experiments  that  document  the  above  flow 
geometry.  Therefore  the  purpose  of  the  present  work  is  to  study  the  behaviour  of  a 
flat  plate  turbulent  boundary  layer  under  the  influence  of  a  circular  cylinder  wake 
meeting  the  leading  edge  at  different  angles.  This  work  also  has  some  relevance 
to  the  influence  of  free  stream  turbulence  on  turbulent  boundary  layers,  with  the 
free  stream  turbulence  being  generated  by  the  wake  of  the  cylinder  rather  than  by 
screens  or  grids  as  in  earlier  work  (reviewed  by  t.g.  Hancock  II  Bradshaw  1983). 

Apparatus 

The  experiments  were  conducted  in  the  910  mm  x  910  mm  low-speed  wind  tunnel  in 
the  Department  of  Aeronautics.  The  experimental  arrangement  is  shown  in  figure 
1.  The  flat  plate  used  in  the  present  study  is  the  same  as  that  used  by  Hancock  & 
Bradshaw  (1983).  A  ogival-profile  nose  piece  with  sharp  leading  edge  was  used  to 
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avoid  intermittent  leading-edge  separation  due  to  the  random  normal-component 
velocity  fluctuations  in  the  wake.  A  nominal  zero  pressure  gradient  along  the  stream 
was  achieved  using  a  flap  to  locate  the  front  stagnation  point  at  the  sharp  leading 
edge.  A  16  mm  diameter  circular  cylinder  was  mounted  ahead  of  the  plate  spanning 
the  test  section.  The  leading  edge  of  the  flat  plate  was  at  0.8  m  from  the  circular 
cylinder,  giving  a  fetch  of  50  diameters  for  the  wake  to  develop  before  it  reaches 
the  plate.  The  circular  cylinder  could  be  held  with  its  axis  at  an  angle  to  the  span 
of  the  flat  plate  (figure  1).  We  refer  to  this  angle  as  the  ‘wake  angle’  for  brevity 
in  our  discussion.  Experiments  were  performed  for  three  wake  angles,  viz.  0°,10° 
and  20°,  at  a  tunnel  velocity  of  16  m/s.  The  measurements  included  skin  friction 
coefficients,  mean  velocity  profiles  and  those  turbulence  quantities  that  appear  in 
the  transport  equations  for  the  Reynolds  stress  components. 

Results 

Figure  2  shows  the  spanwise  skin  friction  distributions  at  different  streamwise  lo¬ 
cations  for  different  wake  angles.  The  distributions  follow  the  expected  trend  in 
that  the  boundary  layer  on  one  half  of  the  surface  (the  right  half  as  seen  looking 
downstream  in  figure  1)  interacts  with  the  wake,  while  the  other  side  of  the  wake 
does  not  seem  to  affect  the  boundary  layer  since  it  passes  under  the  surface.  The 
minima  of  the  distributions  in  the  case  of  non-zero  wake  angle  occur  off  the  centre 
line.  At  large  distances  from  the  centre  line,  the  distributions  tend  to  asymptote  to 
the  two-dimensional  boundary  layer  values. 

Figure  3  shows  the  mean  velocity  profiles  measured  with  a  1.1  mm  diameter  pitot 
tube.  The  profiles  were  measured  in  two  planes,  one  along  the  centre  line  and  the 
other  -80  mm  from  the  centre  line,  roughly  at  the  position  of  minimum  Cj  (figure 
2).  The  streamwise  distance  of  the  profiles  from  the  leading  edge  in  meters  is  as 
x  =  0.152x  (station  number).  The  flat-plate  profile  at  station  14  in  the  absence  of 
the  cylinder  (diamond  symbols  in  figure  3)  demonstrates  how  the  outer  boundary 
conditions  changes  when  the  wake  is  present.  In  fact,  the  reduction  in  C/,  even 
at  0°  wake  angle,  can  be  attributed  to  the  reduced  mean  velocity  as  seen  by  the 
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undisturbed  boundary  layer  without  invoking  any  turbulence  mechanism  as  such. 
The  off-centreline  profiles  at  small  x  retain  the  velocity  minimum  found  in  the 
undisturbed  cylinder  wake,  but  the  profiles  become  monotonic  at  large  x.  Shear 
stress  profiles  are  shown  in  figure  4.  The  off-centerline  profiles  show  some  remarkable 
differences  for  different  wake  angles  due  qualitatively  to  the  change  of  sign  of  shear 
stress  across  the  undisturbed  cylinder  wake.  These  profiles  also  demonstrate  the 
bodily  shift  of  wake  fluid  by  the  growing  boundary  layer,  as  indicated  by  the  upward 
propagation  of  the  maxima  of  the  wake  profile  almost  at  the  same  rate  as  the 
undisturbed  boundary  layer  thickness  with  downstream  distance.  In  addition  all 
the  shear  stress  profile.,  tend  to  extrapolate  to  nearly  the  same  value  at  the  wall 
suggesting  a  weak  influence  of  the  turbulence  in  the  wake  on  the  skin  friction  drag. 

Conclusions 

The  skin  friction  in  a  turbulent  boundary  layer  on  a  flat  plate  decreases  in  the 
presence  of  the  wake  of  an  upstream  mounted  cylinder  representating  an  axial  tur¬ 
bomachine  blade,  in  contrast  to  the  increase  in  Cj  produced  by  multiple- rod  grids. 
Evidently  the  primary  effect  is  the  reduction  in  “free  stream”  velocity  seen  by  the 
boundary  layer,  as  a  consequence  of  the  velocity  defect  in  the  cylinder  wake:  in¬ 
teraction  between  the  turbulence  fields  of  the  wake  and  the  boundary  layer  are 
less  important.  The  details  of  the  flow,  particularly  in  the  present  case  where  the 
cylinder  axis  is  inclined  to  the  plane  of  the  boundary  layer,  are  of  considerable 
interest. 
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EFFECTS  OF  FREE-STREAM  TURBULENCE 
ON  TURBULENT  BOUNDARY  LAYERS 
WITH  CONVECTIVE  HEAT  TRANSFER 

V.  Baskaran#,  O.  E.  Abdellatiff  ii  P.  Bradshaw) 
Dept  of  Aeronautics, 

Imperial  College  of  Science  &  Technology, 

London  Univ. 


ABSTRACT 

Measurements  of  surface  heat  transfer  and  mean 
velocity  have  been  made  in  a  two-dimensional  flat- 
plate  turbulent  boundary  layer,  in  the  presence  of 
grid-generated  free-stream  turbulence.  The  data  show 
that,  as  in  the  case  of  skin-friction  coefficient,  the 
Stanton  number  increases  with  decreasing  turbulence 
length  scale  at  given  turbulence  intensity.  However, 
the  decrease  is  much  more  rapid  than  in  the  skin- 
friction  correlation  suggested  by  Hancock  &  Bradshaw 
(1983).  An  alternative  correlation  is  tentatively  sug¬ 
gested.  Some  comments  are  made  on  previous  data 
and  on  the  use  of  arrays  of  rods  to  simulate  turboma¬ 
chine  blade  wakes. 

NOMENCLATURE 

Cj  -  local  skin  friction  coefficient,  Tw/(l/2pU?) 

Ci  -  constant  in  equation  (3) 

D  -  constant  in  equation  (3) 

Le  -  dissipation  length  scale  (see  equation  2) 

M  -  mesh  size  of  grid 
n  i  -  constant  in  equation  (3) 

R  -  Reynolds  analogy  factor,  St/(Cj /2) 

Ro  -  momentum-thickness  Reynolds  number,  Ue0j v 
St  -  Stanton  number  (see  equation  1) 

Tw  -  wall  temperature 

T„  -  free-stream  temperature 

Ue  -  free-stream  mean  velocity 

u«  -  longitudinal  fluctuating  component  of  free-stream 

#  Aeronautical  Research  Labs.,  DSTO,  Salisbury, 
Australia 

t  Civil  Engg  Dept,  City  Univ.,  London  EC1 
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turbulence 

xg  -  streamwise  distance  from  grid 
x  -  distance  from  plate  leading  edge 
y  -  normal  distance  from  plate  surface. 

AC/  -  incremental  skin  friction  coefficient 

S  -  boundary  layer  thickness,  y  —  6  where  UjUe  = 

0.995 

A  St  -  incremental  Stanton  number 

Suffix  o  -  no  grid,  or  negligible  free-stream  turbulence 

INTRODUCTION 

The  influence  of  free-stream  turbulence  (hereafter 
FST)  on  boundary  layers  is  particularly  relevant  to  the 
flow  over  axial-flow  turbomachine  blades,  which  suffer 
from  the  turbulence  and  periodic  unsteadiness  caused 
by  the  wakes  of  upstream  blades,  leading,  for  instance, 
to  increased  heat  transfer  to  turbine  blades.  It  is  now 
generally  recognised  that  the  effects  depend  not  only 
on  the  intensity  of  FST  but  also  on  its  spectral  distri¬ 
bution  -  most  crudely  represented  by  a  single  length 
scale.  The  present  study  is  the  start  of  a  programme 
on  simulation  of  realistic  turbomachinery  turbulence, 
now  being  continued  at  Stanford. 

The  simplest  form  of  FST,  which  has  attracted  at¬ 
tention  of  many  investigators,  is  that  behind  a  grid, 
since  it  is  easy  to  produce  and  is  close  to  classical  ho¬ 
mogeneous  isotropic  turbulence.  It  was  used  at  the 
start  of  the  present  study  to  check  out  the  techniques 
and  to  improve  knowledge  of  the  effect  of  FST  length 
scale  on  heat  transfer,  of  which  little  is  known  even 
for  the  simplest  case  of  a  flat  plate. 

The  next  stage  of  the  work,  to  be  reported  more 
fully  elsewhere,  was  an  investigation  of  FST  produced 
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by  the  wake  of  a  rod  parallel  to  the  test  plate.  Sev¬ 
eral  investigators  have  used  rotating  or  translating  ar¬ 
rays  of  rods  to  simulate  turbomachine  blade  rows.  We 
started  by  studying  the  effect  of  a  single,  stationary 
rod,  and  found  that  it  decreased  the  skin  friction  on 
the  flat  plate  behind  it:  quite  simply,  the  effect  of  re¬ 
duced  mean  velocity  in  the  rod  wake  was  larger  than 
that  of  the  wake  turbulence!  Because  skin  friction 
varies  roughly  as  (free-stream  velocity)2,  time-average 
skin  friction  for  a  given  mass-flow  rate  is  increased  by 
a  translating  array  of  rods,  simply  because  of  the  in¬ 
creased  free-stream  velocity  between  the  rod  wakes: 
heat  transfer,  which  is  to  a  first  approximation  di¬ 
rectly  proportional  to  free-stream  velocity,  would  be 
comparatively  little  affected.  This  fact  about  rod  ar¬ 
rays  does  not  seem  to  have  been  noticed  by  previous 
users.  Work  on  more  meaningful  simulation  of  tur¬ 
bomachine  turbulence  -  as  distinct  from  periodic  un¬ 
steadiness  -  is  in  progress. 

Several  attempts  have  been  made  to  correlate  chan¬ 
ges  in  skin  friction  coefficient  Cj  and  Stanton  number 
St  with  FST  produced  by  grids.  Simonich  Sc  Brad¬ 
shaw  (1978)  suggested  that  the  increase  in  C /  or  St 
is  directly  proportional  to  the  rms  longitudinal  tur¬ 
bulence  intensity,  with  a  five  percent  increase  in  St 
for  each  one  percent  of  rms  turbulence  (compared  to 
about  3  percent  for  Cj).  They  ignored  any  effects  of 
turbulence  length  scale.  Later,  Hancock  Sc  Bradshaw 
(1983:  herafter  cited  as  HB)  measured  the  effects  of 
grid  turbulence  with  a  wide  range  of  length  scales,  and 
put  forward  a  more  general  form  of  correlation  for  the 
fractional  change  in  skin-friction  coefficient  with  both 
velocity  and  length  scales,  as  a  function  of  the  empir¬ 
ical  variable,  (\fu?JUc)l(Lcl6  +  2).  The  correlation 

i3  strongly  non-linear  in  rms  intensity  ,  and  the 
effects  of  length  scale  are  large,  at  least  in  the  range 
of  the  experiments,  1.5  <  LJ6  <  6.  Blair  (1983)  and 
Castro  (1984)  have  extended  the  correlation  to  low 
Reynolds  numHrs,  where  viscous  effects  on  the  outer 
layer  interact  with  those  of  FST. 

Blair  found  good  agreement  between  his  Stanton- 
number  measurements  in  grid  turbulence  and  the  HB 
correlation  with  the  ordinate  multiplied  by  1.3,  for 
fractional  increments  in  51  up  to  about  0.2.  For  higher 
turbulence  levels,  the  HB  correlation  predicts  that  the 
increase  in  Cy  becomes  slower,  but  Maciejewski  Si 


Moffat  (1988)  and  MacMullin,  Elrod  Sc  Rivir  (1988) 
have  found  that  St  increases  far  more  than  C/  at 
high  levels  of  non-isotropic  FST,  with  roughly  a  fac¬ 
tor  of  four  increase  in  St  at  40  percent  turbulence,  so 
that  Blair’s  conclusion  is  evidently  not  generally  valid. 
(Maciejewski  Sc  Moffat’s  data  do  not  collapse  on  any 
combination  of  length  scale  and  intensity.)  In  1953 
Sugawara,  Sato,  Komatsu  Si  Osaka  (1988)  found  that, 
in  grid  turbulence,  the  factor  of  increase  of  St  tended 
to  an  almost  constant  value  of  about  1.55  at  rms  in¬ 
tensities  of  7  or  8  percent  (the  limit  of  their  data), 
while  Slaniiauskas  Si  PediSius  (1978)  found  an  almost 
constant  value  of  1.4  for  10  to  17  percent  turbulence. 
The  latter  authors  found  that  the  factor  of  increase  of 
St  was  about  1.7  times  the  factor  of  increase  in  Cj  for 
all  FST  values,  but  this  cannot  be  directly  compared 
with  Blair’s  1.3  because  §lan£iauskas  Sc  PediSius  ig¬ 
nored  length-scale  effects.  The  situation  is  confused, 
especially  as  regards  length-scale  effects. 

We  have  repeated  the  experiments  of  Simonich  Sc 
Bradshaw  (1978)  over  about  the  same  range  of  inten¬ 
sities  as  used  by  Blair,  but  over  the  wider  range  of 
length  scales  used  by  HB.  Our  C{  distributions  agree 
with  the  HB  correlation  but  the  Stanton-number  dis¬ 
tributions  do  not,  even  if  the  constant  of  proportion¬ 
ality  is  adjusted:  the  experiment  was  cut  short  by  the 
authors’  departure  from  Imperial  College  but  the  data 
are  offered  as  a  partial  insight  into  length-scale  effects 
over  a  wider  range  than  previous  work.  A  tentative 
modification  to  the  C j  correlation  is  suggested  to  cor¬ 
relate  St  . 

APPARATUS  AND  TECHNIQUES 

The  measurements  were  made  in  the  910  mm  x 
910  mm  closed-circuit  wind  tunnel  used  by  Simonich 
Sc  Bradshaw  (1978)  and  Hancock  Si  Bradshaw  (1983). 
The  flat  plate  (figure  1)  is  also  that  used  by  Simonich 
Si  Bradshaw.  An  ogival  leading  edge  was  used  to  avoid 
intermittent  leading-edge  separation  due  to  normal- 
component  velocity  fluctuations.  Nominally  zero  pres- 
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sure  gradient  along  the  rest  of  the  plate  was  achieved 
by  using  a  trailing-edge  flap.  A  1.2  mm  trip  wire  was 
attached  80mm  from  the  leading  edge.  The  biplane 
grids  used  by  HB,  with  square  bars  at  pitches  of  38 
mm,  76  mm  and  152  mm  were  inserted  at  the  start 
of  the  working  section.  The  leading  edge  of  the  flat 
plate  was  2.09  m  from  the  grid.  All  measurements 
were  made  at  a  free-stream  velocity  of  21  m/s,  as  in 
the  earlier  studies  at  Imperial  College.  At  this  speed 
the  turbulence  level  in  the  empty  working  section  was 
found  to  be  approximately  0.03  percent. 

The  resistive  element  used  to  heat  the  plate  con¬ 
sisted  of  30  stainless  steel  strips,  0.05  mm  thick  and  25 
mm  wide,  connected  in  series  and  sandwiched  between 
two  6  mm  thick  plywood  boards  2.4  m  long,  so  that 
when  steady  conditions  were  reached  the  heat  transfer 
from  each  surface  of  the  plate  was  half  the  electrical 
power  input  and  direct  measurement  was  not  needed. 
The  cold  resistance  of  the  element  was  typically  8fl.  In 
the  present  study,  a  power  input  of  1.3kW  (13A  and 
100V)  was  used:  Simonich  and  Bradshaw  used  only 
450W.  The  mean  temperature  difference  between  the 
surface  and  the  free  stream  was  in  the  range  of  4-7°  C, 
and  was  measured  using  Cromel-Alumel  thermocou¬ 
ples  imbedded  in  the  surface  of  the  plate,  with  the 
cold  junctions  mounted  in  the  free  stream  on  the  tops 
of  Perspex  masts.  The  thermocouples  were  calibrated 
using  ice  blocks  and  hair  driers,  and  the  constant  of 
calibration  was  found  to  be  typically  40 \iV  !aC. 


The  Stanton  number  was  deduced  as 

st  _  _ 

pc„(Tw  -  T.)U. 


(1) 


where  W  is  the  electrical  power  input  and  A  is  the  total 
area  of  the  plate.  Radiation  losses  were  estimated  to 
be  about  10  to  15  percent  of  the  convective  transfer. 

Mean  velocity  profiles  were  measured  on-line  to  a 
BBC  Model  B  microcomputer,  using  a  1.1  mm  pitot 
tube.  Mean  temperature  profiles  were  measured  with 
a  thermocouple  probe,  which  proved  to  suffer  from  ex¬ 
cessive  wall  effects,  so  the  results  are  not  given  here. 
Local  C/  was  deduced  by  fitting  the  mean  velocity 
profiles  to  the  logarithmic  law  of  the  wall  with  con¬ 
stants  0.41  and  5.2,  as  the  law  of  the  wall  was  found 
to  be  unaffected  by  FST  in  the  range  of  intensities 
used  by  HB. 


The  longitudinal  component  of  the  turbulent  in¬ 
tensity  behind  various  grids  was  measured  along  the 
centre-line  of  the  working  section  in  the  absence  of  the 
plate,  using  a  normal  hot-wire  and  a  Melbourne  Uni¬ 
versity  type  constant-temperature  anemometer.  The 
length  scale  of  the  FST  was  derived  from  the  decay  of 
the  longitudinal  component  of  the  turbulent  intensity, 
using  Taylor’s  hypothesis  for  frozen  turbulence,  as 


Lc 


M3/2 

Vc(duVdx) 


(2) 


If  the  turbulence  were  isotropic,  with  turbulent  ki¬ 
netic  energy  k  equal  to  3u*/2,  this  length  scale  would 
be  \/273  times  the  true  dissipation  length  scale  k3^2/t. 

RESULTS 

Figure  2  shows  the  variation  of  the  longitudinal 
turbulent  intensity  for  each  grid,  with  the  limited  data 
of  Simonich  &i  Bradshaw  (1978).  Note  that  the  max¬ 
imum  turbulence  level  that  the  plate  encounters  is 


t'|Cur«  1.  Str.airiwU*  distribution*  of  lli*  longitudinal  turbulent 
intensity  behind  gride. 


about  8  percent.  According  to  Comte- Bellot  Corrsin 
(1966),  the  decay  of  the  longitudinal  turbulent  inten¬ 
sity  can  be  expressed  in  the  form 


where  C|,  D  and  r»i  are  constants  for  a  given  grid 
of  mesh  size,  A/.  Best-fit  values  for  the  present  data 
are  given  in  Table  1. 
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/ 


A/(mm) 

38 

76 

152 

Cf 

0.0659 

0.0476 

0.07 

D 

4.0 

4.81 

3.63 

1.25 

1.25 

1.25 

Table  1. 


The  streamwise  distributions  of  the  length  scale  Le 
derived  from  equation  2  are  shown  in  figure  3,  non- 
dimensionalised  using  the  local  boundary  layer  thick¬ 
ness.  The  length  scale  varies  from  1.5  to  6  times  the 
local  boundary  layer  thickness,  as  in  the  HB  experi¬ 
ments.  Note  that  Le  increases  with  x,  but  more  slowly 
than  6. 
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Figure  4.  Straainwiee  illalrlbtillona  of  Stanton  number  In  tlie  heated 
flat  plate  turbulent  boundary  layer. 


The  streamwise  distributions  of  St  for  the  different 
grids  are  shown  in  figure  4,  which  includes  the  distri¬ 
butions  obtained  by  Simonich  le  Bradshaw  and  in  the 
present  study  without  a  grid.  The  initial  difference 
in  the  St  distribution  between  the  two  cases  is  due  to 
the  adverse  pressure  gradient  aft  of  the  leaaing  edge, 
which  was  stronger  in  Simonich  le  Bradshaw’s  work. 


The  Cf  and  St  distributions  are  plotted  against 
the  momentum  thickness  Reynolds  number  in  figures  5 
and  6,  as  the  effects  of  FST  are  best  evaluated  at  fixed 
R#:  nearly  all  the  results  lie  above  the  upper  limit  of 
significant  interaction  between  the  effects  of  FST  and 
of  low  Reynolds  number  (R#  «  2000  according  to  Cas¬ 
tro).  The  present  Cf  data  corresponding  to  negligible 
FST  (no  grid)  are  in  good  agreement  with  the  general 
Cf-R*  correlation  proposed  by  Coles  (1962)  and  with 
the  earlier  Imperial  College  data.  The  no-grid  Stanton 
lumber  at  R#  =  5000  is  about  4  percent  higher  than 
the  correlation  of  Kays  le  Crawford  (1980). 


Figure  5.  Skin  friction  roellicicuL  vcrmn  momentum  thickness 
Reynolds  mimtier. 
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Figure  7  shows  the  Cf  correlation  of  HB  with  the 
results  of  the  present  and  previous  Imperial  College 
studies.  The  agreement  is  generally  within  ±  2  per¬ 
cent  of  Cf.  The  fractional  changes  in  St  are  plot¬ 
ted  against  the  empirical  parameter  of  HB  in  figure  8, 
together  with  the  limited  data  of  Simonich  le  Brad¬ 
shaw.  Our  measured  values  of  St/St0 ,  as  well  as 
those  of  Simonich  le  Bradshaw,  lie  well  above  the  C/ 
correlation  curve  of  HB.  It  follows  that  the  Reynolds 
analogy  factor  Stj(Cf/ 2)  is  affected  by  FST.  (The 
Reynolds  analogy  factor  is  plotted  against  R*  in  figure 
9  as  an  illustration  of  this  but  the  grid-to-grid  changes 
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Flgur*  8.  Fractional  change*  in  Stanton  number 
ver.u.  ^ • 


in  these  axes  are  not  expected  to  be  simple.)  Exact 
analogy  (St/(C// 2)  =  1)  requires  molecular  and  tur¬ 
bulent  Prandtl  numbers  to  be  unity,  but  in  the  outer 
layer,  where  the  main  effects  of  FST  are  felt,  the  tur¬ 


bulent  Prandtl  number  in  the  absence  of  FST  is  about 
0.7.  With  no  grid,  the  Reynolds  analogy  factor  R  at 
Re  >  3000  is  about  1.22,  but  it  increases  to  roughly 
1.3  at  the  higher  turbulence  levels.  At  lower  Reynolds 
number  the  change  in  St / (C/  /2)  is  larger  but  the  frac¬ 
tional  change  in  St/(C// 2)  —  1  is  roughly  the  same. 

Blair  correlated  his  data  as  R  =  1.18  + 

for  all  length  scales,  a  smaller  variation  than  in  the 

present  data.  It  can  be  shown  easily  that 

ASt/St0  =  (R/R0)[l  +  (AC,/CJa)\  -  1  (4) 

where  R  and  Ra  are  the  Reynolds  analogy  factors 
corresponding  to  finite  FST  level  and  negligible  FST: 
correlations  of  analogy  factor  rather  than  St  may  be 
convenient  for  designers,  but  at  present  it  seems  most 
straightforward  to  modify  the  HB  correlation  for  C / 
and  apply  it  to  St. 
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If  the  data  in  figure  8  are  reliable  they  imply  that 
HB’s  empirical  parameter  for  Cj  underestimates  the 
effect  of  FST  length  scale  on  St.  The  denominator 
would  have  to  be  changed  to  about  (LJ6  -  1)  to  bring 
the  152  mm  grid  results  close  to  the  others,  and  this 
is  clearly  unacceptable  as  a  general  correlation  since 
Le/S  =  1  is  in  the  range  of  practical  interest.  A  fairly 
good  collapse  is  obtained  by  changing  the  denominator 
to  (Le/6  +  2)3,  but  of  course  the  correlation  function 
is  no  longer  a  factored  form  of  Hancock  &  Bradshaw’s. 

CONCLUSIONS 

Skin  friction  coefficients  and  Stanton  numbers  were 
measured  in  a  fiat  plate  boundary  layer  beneath  grid- 
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generated  free  stream  turbulence,  with  a  maximum 
rms  intensity  of  8  percent  of  the  free  stream  veloc¬ 
ity  and  a  maximum  length  scale  of  6  times  the  local 
boundary  layer  thickness.  The  increase  in  skin  friction 
coefficient  agrees  well  with  the  correlation  curve  based 

on  the  variable,  \J^/Ue/(Le/6- t-2),  proposed  by  Han¬ 
cock  ( k  Bradshaw  (1983).  Stanton  number  increases 
do  not  seem  to  follow  the  proposed  skin  friction  cor- 
,  relation  or  any  simple  multiple  of  it,  indicating  that 
the  Reynolds  analogy  factor  is  changing  significantly  - 
more  than  in  the  work  of  Blair  (1983).  A  change  in  the 
above  empirical  variable  to  improve  Stanton  number 
predictions  is  tentatively  suggested. 

The  wake  of  a  spanwise  rod  mounted  ahead  of  the 
test  plate  to  simulate  turbomachine  blade  wakes  pro¬ 
duced  a  reduction  in  skin  friction,  calling  into  question 
the  interpretation  of  published  results  from  moving 
rod  arrays. 
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